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ABSTRACT 

We report the results of an X-ray observing campaign on the massive, evolved star rj Car, 
concentrating on the 2003 X-ray minimum as seen by the XMM-Newton observatory. These are 
the first spatially-resolved X-ray monitoring observations of the stellar X-ray spectrum during 
the minimum. The hard X-ray emission, believed to be associated with the collision of r\ Car’s 
wind with the wind from a massive companion star, varied strongly in flux on timescales of days, 
but not significantly on timescales of hours. The lowest X-ray flux in the 2 — 10 keV band seen 
by XMM-Newton was only 0.7% of the maximum seen by RXTE just before the X-ray minimum. 
In the latter half of the minimum, the flux increased by a factor of 5 from the lowest observed 
value, indicating that the X-ray minimum has two states. The slope of the X-ray continuum 
above 5 keV did not vary in any observation, which suggests that the electron temperature of 
the hottest plasma associated with the stellar source did not vary significantly at any phase. 
Through the minimum, the absorption to the stellar source increased by a factor of 5 — 10 to N h 
~ 3~4xl0 23 cm” 2 . The thermal Fe XXV emission line showed significant excesses on both the 
red and blue sides of the line outside the minimum and exhibited an extreme red excess during 
the minimum. The Fe fluorescence line at 6.4 keV increased in equivalent width from 100 eV 
outside the minimum to 200 eV during the minimum. The small equivalent widths of the Fe 
fluorescence line suggests small fluorescence yield in the companion’s low-density wind. The lack 
of variation in the plasma temperature is consistent with the eclipse of the X-ray plasma during 
the minimum, perhaps by a clumpy wind from the primary star, although the deformation of 
the Fe XXV profile and the relatively weak fluorescence Fe line intensity during the minimum 
may suggest an intrinsic fading of the X-ray emissivity. The drop in the colliding wind X-ray 
emission revealed the presence of an additional X-ray component which exhibited no variation on 
timescales of weeks to years. This new component has relatively cool temperature (kT ~ 1 keV), 
moderate N& (^5xl0 22 cm” 2 ), large intrinsic luminosity (Lx ^10 34 ergs s” 1 ) and a size <1" 
(2300 AU at 2.3 kpc). This component may be produced by the collision of high speed outflows 
at v ~ 1000— 2000 km s” 1 from rj Car with ambient gas within a few thousand AU from the star. 
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1. Introduction 

77 Car is a violently unstable, extremely lumi- 
nous object and a key tracer of evolution of stars in 
the upper portion of the Hertzsprung-Russell dia- 
gram (Davidson & Humphreys 1997). The star is 
believed to have had an initial mass of >150 M 0 
(Hillier et al. 2001). It is currently in a short, 
poorly understood evolutionary stage, known as 
the Luminous Blue Variable (LBV) phase, which 
is thought to occur near the onset of pulsational 
instabilities. It provides a convenient laboratory 
to study how extremely luminous, massive stars 
evolve and how they shape their environments 
both geometrically and chemically. 

77 Car is best known for an extraordinarily pow- 
erful eruption in 1843 which sent ~12 M 0 (Smith 
et al. 2003b) of its atmosphere into space, creat- 
ing a beautiful bipolar nebulosity called the Ho- 
munculus around the star, 77 Car also had a minor 
eruption in 1890, which produced a small bipolar 
nebula inside the Homunculus nebula (the “Little 
Homunculus”, Ishibashi et al. 2003). The star still 
exhibits a strong mass loss (10“ 4 — 10~ 3 M 0 yr”" 1 
Davidson et al. 1995; Cox et al. 1995; Ishibashi 
et al 1999; Hillier et al. 2001; Pittard & Corcoran 
2002), preferentially in the polar direction (Smith 
et al. 2003a; van Boekel et al. 2003). The ejecta 
and the stellar wind are rich in helium and nitro- 
gen and depleted in oxygen and carbon (Verner 
et al. 2005; Davidson et al. 1984: Hillier et al. 
2001). consistent with CNO processing. Davidson 
et al (2005) and Martin & Koppelman (2004) sug- 
gested from the change of Ha and H fi line profiles 
and the increase of optical brightness that 77 Car 
is changing rapidly at the present time. 

The past decade has witnessed an important 
change in our understanding of the star. This 
began with the recognition that the strength 
of some narrow emission lines, notably He I 
10830 A, vary predictably with a period of 5.52 
years (Damineli 1996). Along with observations 
in infrared (Whitelock et al. 1994; Feast et al. 
2001), mm (Cox et al. 1995), cm (Duncan et al. 
1995), and optical (van Genderen et al. 1999) 
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wavelengths, observations of the X-ray emission 
from the star have played a key role. ROSAT ob- 
servations in 1992 first showed a variation which 
appeared to be correlated with Damineli’s emis- 
sion line variations (Corcoran et al. 1995), and 
subsequently the 2—10 keV lightcurve of the 
star obtained by RXTE (Ishibashi et al. 1999) 
showed in detail the X-ray variation culminating 
in a swift, unstable rise to maximum and steep fall 
to minimum lasting for ~3 months, i.e. the same 
time interval as Damineli’s spectroscopic mini- 
mum. These observations clearly show that the 
star varies in a fundamental way every 5.54 years, 
and strongly suggest that 77 Car is two stars, not 
one. In this binary model, most of the dramatic, 
pan-chromatic changes are now believed to be pro- 
duced by the interaction of the UV flux and wind 
from a companion star with the wind of 77 Car. 
Variations in the X-ray region are produced by 
a wind-wind collision (WWC) of the wind from 
77 Car with the wind of the companion star. A 
current guess at the system parameters describes 
a massive hot companion with M c ~30 ikf 0 plus a 
brighter primary star with M n >80 M 0 in a highly 
eccentric (e >0.9), 5.54 year orbit (Corcoran et al. 
2001a, see also Davidson 1999). 

During the 1997/98 X-ray minimum, 77 Car 
was the target of every available X-ray obser- 
vatory, including RXTE (Ishibashi et al. 1999; 
Corcoran et al. 2001a), ASCA (Corcoran et al. 
2000), and BeppoSAX (Viotti et al. 2002). The 
ASCA satellite detected hard X-ray emission from 
77 Car (Corcoran et al. 2000), which was char- 
acterized by slightly smaller iVn than the pre- 
minimum state with reduced plasma emission 
measure {EM.). The limited spatial resolution 
of ASCA , however, left the possibility that the 
observed emission was due to contamination by 
unresolved nearby sources. The BeppoSAX satel- 
lite observed 77 Car just after the recovery of the 
X-ray emission (Viotti et al. 2002). The spectrum 
showed strong absorption of Ah~1-5x10 23 cm* 2 
unlike the ASCA spectrum during the minimum, 
which had an Ah ^3x10 22 cm“ 2 . The overall X- 
ray brightness variations observed with RXTE are 
explained well with the colliding wind mechanism 
(Pittard et al. 1998; Pittard & Corcoran 2002). 
However, key properties of the X-ray lightcurve 
were still poorly understood: dramatic changes 
in X-ray flux as the emission increases to maxi- 
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Table 1 

Logs of the XMM-Newton and Chandra Observations 


Abbr. 

Seq. ID 

Observation Start 
Date 

4>orbit 

Exposure 

(ksec) 

Observation Mode 

Filter 

Before the minimum 
XMM 00 0726 a 112580601 

2000 Jul. 26, 06:01 (M2) 

0.470 

••• / /29.6 

PFW/PFW/PSW 

thick 

XMM 00 o72S a 

112580701 

2000 Jul. 28, 00:51 (M2) 

0.471 

/ ••• /7.8 

PFW/PFW/PSW 

thick 

XMMo30125 b 

145740101 

2003 Jan. 25, 13:03 (PN) 

0.922 

4.8/ / ... 

PSW/PFW/PFW 

thick 

XMMo30127A b 

145740201 

2003 Jan. 27, 01:08 (PN) 

0.922 

4.8/ ... / 

PSW/PFW/PFW 

thick 

XMM 03 O X27B b 

145740301 

2003 Jan. 27, 20:42 (PN) 

0.923 

4 . 7 / ••• / ... 

PSW/PFW/PFW 

thick 

XMMo30129 b 

145740401 

2003 Jan. 29, 01:45 (PN) 

0.923 

5 . 8 / ... / ... 

PSW/PFW/PFW 

thick 

XMM 03 0130 b 

145740501 

2003 Jan. 30, 00:00 (PN) 

0.924 

4.8/ ... / ... 

PSW/PFW/PFW 

thick 

XMM 030608 

160160101 

2003 Jun. 08, 13:31 (Ml) 

0.988 

22.1/27.6/ 

PSW /PSW/PFW 

thick 

XMMq 30613 

160160901 

2003 Jun. 13, 23:52 (Ml) 

0.990 

21.9/30.4/ 

PSW/PSW/PFW 

thick 

During the minimum 
CXO 030720 c 200216 

2003 Jul. 20, 01:47 

1.008 

90.3 

ACIS-S w/HETG 


XMMo 30722 C 

145780101 

2003 Jul. 22, 01:51 (M2) 

1.009 

5.8/8.2/S.4 

PSW/PSW/PFW 

thick 

XMMo 30802 

160560101 

2003 Aug, 02, 21:01 (M2) 

1.015 

12.2/17.1/17.5 

PSW/PSW/PFW 

thick 

XMMq 30 S 09 

160560201 

2003 Aug. 09, 01:44 (M2) 

1.018 

S. 7/12.3/12. 6 

PSW/PSW/PFW 

thick- 

XMMo 30 S 1 S 

160560301 

2003 Aug. 18, 15:23 (M2) 

1.023 

13.0/18.2/18.6 

PSW/PSW/PFW 

thick 

CXOo30828 

200237 

2003 Aug. 28, 17:38 

1.028 

18.8 

ACIS-S 


After the minimum 
CXO 030926 200217 

2003 Sep. 26, 22:47 

1.042 

70.1 

ACIS-S w/HETG 



Note. — Abbr.: abbreviation adopted for each observation. Seq. ID: sequence identification number of each obser- 
vation. Observation ' Start Date: start time of the detector in the parentheses. <f>orbit = (JD[observation start] — 
2450799. 792)/2024 (Corcoran 2005). Exposure: exposure time excluding the detector deadtime; for XMM-Newton ob- 
servations the exposure times are given for the PN, MOS1 and MOS2 detectors, respectively; The data sets without 
exposure time were not used due to the severe event pile-up. Observation Mode: PSW: prime small window. PPW: 
prime full window. For the XMM-Newton observations the observations modes are given for the PN, MOS1 and MOS2 
detectors, respectively. Filter: optical blocking filter selected for the XMM-Newton detectors. 

a Data from XMM000726 and XMM 000728 are combined and denoted as XMMqoo 7 - 

b Data from XMMo30i25t XMMo3oi27a, XMMo30127B> XMM 030129 and XMM 030130 are combined and denoted as 
XMMq301 . 

c Spectra from CXO030720 and XMM030722 are simultaneously fit and those results are denoted as CXOXMMq307- 
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mum, the variation of the absorbing material in 
front of the X-ray emitting region near the X-ray 
minimum, the nature of the rapid fall from X-ray 
maximum to X-ray minimum, the excess in Ah 
after recovery. 

The minimum which occurred in mid- 2003 was 
among the best observed astronomical events of 
all time. A key part of this campaign were the 
detailed X-ray observations obtained with RXTE , 
Chandra , and XMM-Newton. RXTE again pro- 
vided crucial monitoring of the daily changes in 
X-ray flux, while Chandra provided previously- 
unobtainable monitoring of the X-ray emission 
line dynamics at key phases of the binary period. 
XMM-Newton obtained critical measurements of 
the spectrum of the source during the low-flux 
state when observations with RXTE or Chandra 
are difficult. 

This paper describes the overall change of 
the X-ray spectrum of 77 Car around the 2003 
minimum as measured by XMM-Newton (Jansen 
et al. 2001), supplemented where necessary with 
key observations with Chandra (Weisskopf et al. 
2002). We try to understand the X-ray emission 
and absorption mechanism and the cause of the 
X-ray minimum, from a comparison of derived 
spectral properties such as kT ) Ah and emission 
line strengths in these observations and earlier 
ROSAT and A SC A observations during the 1992 
and 1997/98 minima. 

The paper is comprised of the following sec- 
tions. Section 2 describes the observations and 
method of the data reduction, and includes a de- 
scription of problems of unresolved emission com- 
ponents around the central source in the extracted 
XMM-Newton source events. Section 3 describes 
imaging and timing analyses of the XMM-Newton 
data. Section 4 describes the analysis of the X-ray 
spectra during the 2003 minimum. Section 5 com- 
pares the X-ray spectra during the 2003 minimum 
to the previous minima in 1998 and 1992. We dis- 
cuss our results in section 6, and in section 7 we 
summarize our conclusions. 

2. Observations & Data Reduction 

In 2003, XMM-Newton observed 77 Car a total of 
11 times, 5 times in January 2003 (about 6 months 
before the onset of the 2003 minimum), twice in 
June 2003 (weeks before the onset of the 2003 min- 


1381/92 Dccl Jaril Febl Marl Aprl Mayl Jual Jull Aujjl Scipl 

1997/88 Jtill Augl Scjpl Octl Novi Peel Janl Fobl Marl Aprl 


2003 Jaal Febl Marl A|>rl Mayl Junl Jull Augl Scjpl Octl 



Fig. 1. — RXTE light curves (solid lines, black : 
2003, red: 1997/98) around the X-ray mini- 

mum with the timings of XMM-Newton , Chandra , 
ASCA, BeppoSAX , and ROSAT pointed observa- 
tions (black: 2003, red: 1997/98, green: 1992). 
Multiple ROSAT observations between 1.00 to 
1.01 were summed together for ROSAT92oei2- 

imum, near the X-ray maximum) and four times 
during the 2003 minimum in July and August 
2003 prior to the recovery, which started near 2003 
September 3. Figure 1 displays these observing 
times on the RXTE lightcurve (Corcoran 2005) for 
context. We also re- analyzed two XMM-Newton 
observations from 2000, whose results were previ- 
ously published by Leutenegger et al. (2003). We 
thus used a total of 13 XMM-Newton observations 
in our analysis, as listed in Table 1. Throughout 
this paper, individual XMM-Newton observations 
are designated XMM, subscripted with the year, 
month and day of the observation. To improve 
statistics, we combined data from two observations 
in 2000 and combined data from five observations 
in 2003 January because these observations were 
close together in time, and did not show any sig- 
nificant variabilty among them. 

XMM-Newton is composed of three nested 
Wolter I- type X-ray telescopes (Aschenbach et al. 
2000) with the European Photon Imaging Camera 
(EPIC) CCD detectors (PN, MOS1 and MOS2) 
in their focal planes (Struder et al. 2001; Turner 
et al. 2001). In all the observations, 77 Car w*as 
on-axis but the observations were obtained with 
different satellite roll angles. The EPIC instru- 
mental modes are listed in Table 1. The XMM- 
Newton spectra suffered significant photon pileup 
in the data taken in prime full window (PFW) 
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mode outside the 2003 minimum when the star 
was bright in the 2-10 keV band. We did not use 
these data sets in our analysis. XMM-Newton is 
also equipped with the Reflection Grating Spec- 
trometer (RGS; den Herder et al. 2001), but the 
RGS has limited sensitivity above ^1 keV, where 
emission from the WWC region dominates. We 
therefore did not use the RGS data. 

The analysis of the XMM-Newton data was per- 
formed with version 5.4.1 of the SAS 1 software 
package and version 5.2 of the HEASoft 2 analy- 
sis package. The Observation Data Files (ODF) 
data were processed using the SAS scripts “em- 
chain” and “epchain.” We removed events close to 
hot pixels or outside the field of view, and selected 
events with pattern ^4 for EPIC PN spectral anal- 
ysis and pattern ^12 for timing analysis of the 
EPIC PN data and spectral and timing analysis of 
the EPIC MOS data. Most of these observations 
luckily avoided high background periods and the 
instrumental background is negligible in any ob- 
servations within the 0.3—10 keV band. We did 
not need to reject any high background periods 
because the strong emission from rj Car greatly 
surpassed the observed background levels. 

Due to the limited spatial resolution of XMM- 
Newton, the extracted 77 Car data inevitably in- 
clude emission from a number of sources (see Fig- 
ure 2). The extracted spectrum is usually domi- 
nated by hard X-rays from the wind- wind collision 
in r\ Car, which varies from observation to observa- 
tion (Corcoran et al. 2000; Viotti et al. 2002). The 
Homunculus Nebula emits weak diffuse reflected 
X-rays from the central source, which can be spa- 
tially resolved only with Chandra and only during 
the minimum (Corcoran et al. 2004). The outer 
ejecta extending about 1' beyond the Homuncu- 
lus Nebula emits diffuse non- variable X-rays below 
~1 keV with a ring- or shell-like morphology (Se- 
ward et al. 2001; Weis et al. 2004). To correct the 
XMM-Newton spectra for emission from the outer 
ejecta and the Homunculus Nebula, we used three 
Chandra data sets from the 2003 X-ray observ- 
ing campaign, two observations during the 2003 
minimum and an observation just after the 2003 
minimum. These Chandra observations spatially 
resolve the outer ejecta emission from the emission 


x http:/ /xmm. vilspa.esa.es/external/xmm_sw_cal/sas_frame.shtml 
2 http://heasarc.gsfc.nasa.gov/ docs / software/lheasoft / 


of the stellar source and so allow us to determine 
the spectrum of the outer emission, which we can 
then use to help model the XMM-Newton spectra. 

The log of Chandra observations and the tim- 
ing of the observations are also given in Table 1 
and Figure 1 . Chandra observations are desig- 
nated CXO, subscripted with the year, month 
and day of the observation, similar to our des- 
ignation of the XMM-Newton observations. The 
Chandra observations were obtained using the Ad- 
vanced CCD Imaging Spectrometer detector us- 
ing the Spectrometer array (ACIS-S) either with 
or without the high energy transmission grating 
(HETG). For the grating data, Oth-order photon 
events were used for CXO030720 (during the X- 
ray minimum, when the source was too faint at 
1st order), and lst-order photon events were used 
for CXOo30926 after the recovery since the central 
source had severe pile-up in the Oth-order data. 
The observation CXOo3082$, taken during the min- 
imum with the ACIS-S imaging array with no grat- 
ing, also suffered mild pile-up (~ 4%). The anal- 
ysis of the Chandra data was performed with the 
CIAO 3 software package, version 2.3 (CALDB ver- 
sion 2.22), version 3.0 (CALDB version 2.23) and 
version 3.0.2 (CALDB version 2.26) for sequences 
200216, 200237, and 200217 respectively. We also 
used version 5.2 of the HEAsoft package, and fol- 
lowed the recommendations of the analysis science 
thread 4 . A full analysis of the Chandra spectra 
will be published separately (Corcoran et al. 2007, 
in preparation). 

3. X-ray Images and Time Variability 

All XMM-Newton images similarly showed 
emission from the central source and the outer 
ejecta (Figure 2). Because these XMM-Newton 
images cannot completely separate emission from 
the outer ejecta and central source nor resolve 
the faint X-ray Homunculus Nebula (Corcoran 
et al. 2004), we used an ellipse of 50" x 37."5 as 
the source region (called “XMM SRC” in Figure 
2), which includes all these emission components 
but excludes nearby X-ray luminous stars (such 
as HDE 303308) more than ~ V distant. For 
some small window mode MOS observations, the 
south-east tip of the ellipse fell outside the small 
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Fig. 2. — True color images of r\ Car. left: Chandra HETG 0th~order data from CXO 030720 (red: 0.2—1 keV, 
green: 1-3 keV, blue: 3—10 keV); middle : XMM-Newton image accumulated from all the available MOS 
data taken in prime full window mode (red: 0.2-1 keV, green: 1-3 keV, blue: 3-9 keV); right: ASCA GIS 
data from ASCA 971224 (red: 0.8-1 keV, green: 1-3 keV, blue: 3-10 keV). The color in the ASCA image 
has been chosen to highlight the diffuse emission. Solid lines show source (SRC) or background (BGD) 
regions (Designations: OE: Outer Ejecta, HN: Homunculus Nebula, CS: Central Source). The Chandra and 
ASCA images were smoothed with a Gaussian with a — 1.5 pixel, while the XMM-Newton image was not 
smoothed. 


window. We determined background using regions 
from apparently source-free regions on the same 
CCD chip. 

Table 2 shows detected photon count rates in 
the soft (0.3—1 keV), medium (1-4 keV), and hard 
(4-10 keV) bands in each detector for each obser- 
vation. The soft band mainly includes emission 
from the outer ejecta, while the medium band in- 
cludes emission from an additional source which 
we call the central constant emission component 
(CCE, see §4.1.2), and the medium band out- 
side the minimum and the hard band includes the 
WWC emission (see §4.2). In the soft band, the 
MOS2 count rates are 10% lower in XMM 0007 than 
during the 2003 minimum because of the use of 
a smaller MOS2 source region for the XMM 0007 
data. The soft-band PN count rates are slightly 
higher outside minimum because the spectral re- 
sponse of the PN allows photons from the bright 
hard source to contaminate the soft band. For 
example, we estimate the soft band contamina- 
tion due to the hard photons in XMM 0007 to be 
~0.08 counts s'" 1 , which makes the corrected PN 
count rates in the soft band 1.06 counts s”" 1 in this 
observation. This means that the emission from 
the outer ejecta was constant to <-2% to better 
than 90% confidence over a half-year interval. On 


the other hand, both the medium and hard band 
count rates decreased dramatically to ~l-5% of 
the maximum brightness seen by XMM-Newton . 

We constructed light curves with 500 s time 
bins, combining all the available MOS and PN 
data. In Figure 3, we show three sample light 
curves, representing variations before the min- 
imum (XMM 030613 ) and during the minimum 
(XMMo 30722 , XMM 0308 I 8 ) • In general, these light 
curves showed no significant variation in most 
energy bands at >90% confidence (Table 2, Fig- 
ure 3). The exceptions are the medium band light 
curves before the 2003 minimum. The medium 
band light curve from the XMM 0306 O 8 observation 
showed an “excess” of ~0.2 counts s” 1 (~4 %) 
near the middle of the observation, while the 
medium band light curve from the XMM 030613 
observation showed a slight linear increase by 
~ 0.54^22 c °unts s"” 1 day"" 1 (l~cr) or ^10% 
day*" 1 . The linear model is still not accepted 
above the 90% confidence level (x 2 /d.o.f. — 

1.59), due to small fluctuation on timescales of 
<500 s. We note that the hard band light curve in 
XMM 030613 is also better reproduced by a linearly- 
increasing model with a slope of ~ 0.78lto;26 
counts s'" 1 day”" 1 , or ^7.5% day”" 1 , with x 2 /d.o.f. 
= 0.98. This increase is roughly consistent with 
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Fig. 3.-— Background subtracted EPIC PN + MOS light curves from XMM 030613 (top), XMM 030722 ( bottom 
left), and XMM 0308 I 8 ( bottom right). Top and bottom graphs in each panel show the medium (1-4 keV) and 
hard (4-10 keV) bands, respectively. A solid line in each panel shows the best-fit constant model. 
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Table 2 

Time Variability 


Observation Net Count Rate Constant Model Fit 



Soft[0.3— 1 keV] 

PN MOS1 MOS2 
(counts s -1 ) 

Med [1—4 keV] 

PN MOSl MOS2 

(counts s“ x ) 

Hard [4— 10 keV] 

PN MOSl MOS2 

(counts s” 1 ) 

Soft 

X 2 /d.o.f. 
Med Hard 

d.o.f 

XMM 0007 



0.255 a 



1.25 a 



0.91 a 

0.97 

1.12 

0.82 

76 

XMMosox 

1.15 



7.72 



7.18 



0.98 

1.14 

1.05 

75 

XMMososos 

1.10 

0.247 a 


4.58 

1.70 a 


6.77 

1.89 a 


0.99 

1.45 

0.95 

55 

XMMo30623 

1.09 

0.242 a 


4.99 

1.84 a 


10.49 

2.80 a 


1.02 

1.84 

1.29 

62 

XMM 030722 

1.03 

0.223 a 

0.285 

0.456 

0.170 a 

0.188 

0.183 

0.045 a 

0.049 

0.81 

0.85 

1.31 

16 

XMM 030802 

1.07 

0.222 a 

0.288 

0.464 

0.161 a 

0.190 

0.517 

0.137 a 

0.124 

0.97 

1.21 

0.75 

34 

XMM 030 S 09 

1.06 

0.226 a 

0.280 

0.476 

0.178 a 

0.188 

0.554 

0.159 a 

0.151 

1.13 

0.79 

1.09 

24 

XMM03081S 

1.05 

0.22l a 

0.288 

0.514 

0.169 a 

0.202 

0.648 

0.179 a 

0,174 

1.34 

0.96 

0.92 

37 


^Source regions are slightly smaller for these Prime Small Window observations. 


the RX TE light curve 1 5% day 1 , see F igure 1) . 

4. X-ray Spectra 

For each XMM-Newton observation, we pro- 
duced EPIC PN and MOS spectra using the same 
source and background regions as in §3. Figure 
4 presents the EPIC PN spectra for all 2003 ob- 
servations. The X-ray spectra are a combination 
of a variable hard WWC component, along with 
non- variable emission from a number of compo- 
nents. Consistent with Table 2, the spectra below 
1 keV are almost constant through the observa- 
tions. Above 1 keV the spectra varied by a factor 
of 2-3 before the 2003 minimum, then decreased 
dramatically by more than an order of magnitude 
after the onset of the. 2003 minimum in June, with 
small recovery by a factor of ~3 at E >4 keV dur- 
ing August. Most spectra clearly showed lines of 
hydrogen- and helium- like ions of Mg, Si, S, Ar, 
Ca, Fe, and Ni, and a fluorescence line from cold 
Fe. 

In the following we examine first the non- 
variable emission, then use this emission to fully 
understand changes in the variable component. In 
our analysis, we adopt abundances relative to solar 
abundances given by Anders & Grevesse (1989) 

4.1, Non- Variable Components 

4 >1.1. Contribution of the Outer Ejecta and Ho- 
munculus Nebula to the X-ray Spectra 

We use Chandra observations to estimate the 
amount of contamination from the outer ejecta 


and the Homunculus Nebula in the XMM-Newton 
spectra. 

We extracted the spectrum of the outer ejecta 
from the Chandra CXO 030720 HETG Oth-order 
data, which has the smallest contamination from 
the central hard source. We extracted the source 
spectrum from an 83" x 47" ellipse centered on the 
star (“OE SRC” in Figure 2), with its major axis 
parallel to the long axis of the extended outer X- 
ray nebula, but excluded a 20" x 15" ellipse (“HN 
SRC” in Figure 2) centered on the central source 
and the Homunculus Nebula. We extracted back- 
ground spectra from a nearby apparently source- 
free region, rj Car is embedded in soft diffuse emis- 
sion from the Carina nebula, but its contribution 
to the 77 Car spectrum is negligible. 

The spectrum from the “OE SRC” region, ex- 
cluding the spectrum from the “HN SRC” re- 
gion, is shown in Figure 5. This spectrum shows 
a strong emission line from nitrogen, which re- 
confirms the N enhancement in the outer ejecta 
found by a number of previous analyses, (opti- 
cal/UV: Davidson et al. 1982, X-rays: Tsuboi 
et al. 1997 [N/O ~5], Leutenegger et al. 2003 [N/O 
>9]). This spectrum also shows a significant ex- 
cess at 1.3 keV, the energy of the helium-like 
Mg line, which seems to require a Mg overabun- 
dance of >3 solar at kT M3.6 keV (though plasma 
with that temperature should emit a stronger 
hydrogen-like Mg line). Leutenegger et al. (2003) 
also found a similarly strong Mg line in the XMM- 
Newton grating spectrum, though they conjec- 
tured that some of the Mg emission might be con- 
tamination from the stellar X-ray source. Our re- 
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0.5 1 2 5 10 

channel energy (keV) 


Fig. 4. — XMM-Newton EPIC PN spectra of r) Car in 2003. The central source [central constant emission 
(CCE) and colliding wind emission (CWE)] component, the Homunculus Nebula (HN) and the outer ejecta 
components (OE) are marked. Identified emission lines are marked by lines at their rest energies. Fe I marks 
the location of the Fe absorption edge. The red, blue, and grey lines show Chandra best-fit models of the 
outer ejecta (see Figure 5), the Homunculus Nebula, and their sum, folded through the EPIC PN effective 
area. 
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suits suggest that the Mg line does in fact orig- 
inate in the outer nebulosity. We modeled the 
extracted, background-corrected source spectrum 
below 3 keV (where the emission from the ex- 
tended source dominates) using a simple absorbed 
optically thin, thermal plasma model (WABS, 
Morrison & McCammon 1983; MEKAL, Mewe 
et al. 1995), including the ACISABS 5 compo- 
nent to compensate for the progressive low en- 
ergy degradation of the quantum efficiency of the 
Chandra ACIS detector. We allowed the abun- 
dances to vary, but constrained all abundances to 
have the same value except for nitrogen, which 
was allowed to vary independently. We added a 
Gaussian with a fixed line centroid at 1.31 keV 
with no intrinsic broadening, to account for the 
helium-like Mg line emission. The model suc- 
cessfully reproduces the spectrum of the outer 
ejecta, yielding a reduced x 2 of 1.09 for 81 dof, 
with the best-fit parameters of kT ~0.58 keV, Ah 
~ 7.2 xlO 20 cm" 2 , with a VMEKAL normalization 
of 7.2 x 10 ~ 4 , elemental abundance of ~85 solar for 
nitrogen, near solar abundance for the others, and 
an excess flux of ~1.5 xlO"* 5 counts cm“ 2 s^ 1 for 
Mg. The plasma temperature, ~ 0.6 keV - mainly 
constrained by continuum above 1 keV - agrees 
well with the highest plasma temperature derived 
from the emission lines in the XMM-Newton grat- 
ing spectrum (Leutenegger et al. 2003). However, 
the Ah is smaller than interstellar absorption to 
rj Car (>3 x 10 21 cm“ 2 , see §2.2 of Leutenegger 
et ah 2003), the nitrogen abundance is overesti- 
mated from earlier measurements, and the best- fit 
model does not reproduce the Fe L and Ne lines 
near 1 keV that are observed with the high reso- 
lution grating spectrum. This is perhaps caused 
by assuming a simplistic IT model for a spectrum 
which may have multiple cooler plasma compo- 
nents. Therefore, this model reproduces CCD- 
resolution spectra (>100 eV at 1 keV) at energies 
above ~0.4 keV but not higher-resolution grating 
spectra. 

The reflected X-ray emission from the Ho- 
munculus Nebula did not vary significantly be- 
tween CXO 030720 and CXO Q 30 S 2 S ? an interval of 
40 days (Corcoran et al. 2004). Since this emission 
represents the accumulated X-rays reflected from 


5 http : / / cxc . harvard . edu / cont- 
soft / software/ ACISABS. 1 . 1 .html 
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channel energy (keV) 

Fig. 5. — Chandra spectrum of the outer ejecta 
from CXO 03 0720 (OE SRC - HN SRC in Fig- 
ure 2). The solid line shows the best-fit absorbed 
IT model and the dashed line depicts the Gaussian 
component for helium- like Mg. 

the entire Homunculus Nebula, any time varia- 
tion of the central source shorter than the average 
light-travel time from the star to the reflecting 
site inside the Homunculus Nebula (~88 days) 
will be smeared out. We assume that the emission 
from the Homunculus Nebula did not vary be- 
tween CXO 030720 and CXO 03 Q 828 ? which includes 
all XMM-Newton observations during the 2003 
minimum (XMM030722 “ XMMo3o$i8)* To show 
the relative magnitude of the contamination from 
this reflected component, in Figure 4 we include a 
model fit of the Chandra spectrum of the reflected 
X-rays from the Homunculus Nebula, which we 
extracted from a 20" x 15" ellipse centered on 
the star (HN SRC in Figure 2), excluding a 2."5 
radius circle (CS SRC). 

Though the X-ray intensity of the Homuncu- 
lus in observations outside the 2003 minimum is 
less certain, at these times the X-ray emission 
from the central source was so bright that the X- 
ray contamination from the Homunculus is <1% 
even if the reflected emission from the Homuncu- 
lus brightened by a factor of ~3. 

4-1-2. Discovery of a “Constant” Component 
Near the Central Source 

Figure 4 depicts the XMM-Newton EPIC PN 
spectra recorded through 2003, along with the 
best-fit spectral models of the outer ejecta, the 
X-ray Homunculus Nebula, and their sum. Ex- 
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Table 3 

Best Fit Parameters op the CCE Component 


Parameters 


IT fit 

kT 

(keV) 

1.05 (0.96-1.13) 

Zt Si) 

(solar) 

0.22 (0.14-0.29) 

Z{ S) 

(solar) 

0.47 (0.39-0.55) 

E.M. 

(10 56 cm~ 3 ) 

7.0 (6.0-8.2) 

Nk 

(10 22 cm~ 2 ) 

4.96 (4.77-5.16) 

Lx & 

(IQ 34 ergs s" 1 ) 

1.3 

X 2 /d.o.f 


1.51 

d.o.f 


1205 


NOTE. — Elemental abundances of the hard com- 
ponent were found to be non-solar, but do not vary 
with observation. 

a Absorption corrected X-ray luminosity between 
0.5-10 keV, assuming the distance of 2.3 kpc. 


cept for the excess below 1 keV, which is produced 
by the poor absolute flux calibration of the EPIC 
PN 6 , the excess over the summed spectrum repre- 
sents emission in a region (CS SRC in Figure 2) 
within 2/ / 5 of the hard X-ray peak. Apparent in 
Figure 4 is an excess over the model in the energy 
range 1 — 3 keV. This emission apparently did not 
vary in time with the harder emission, and had 
the same flux in the Chandra spectra of the central 
region in observations CXO 030720 , CXO 030828 ? and 
CXO 030926 (see also Figure 6), so that the emission 
was securely constant for more than two months 
between July 20 - September 26 2003. We call this 
component the Central Constant Emission compo- 
nent, or “CCE” . 

The CCE shows lines of helium-like S and Si 
ions but no lines of hydrogen-like S and Si ions. 
This indicates that the CCE comes from relatively 
cool plasma. To estimate the emission parameters 
of the CCE, we subtracted contributions of the 
outer ejecta and Homunculus Nebula from both 
the PN and MOS spectra (see the top panel of 
Figure 6 which shows the PN spectra). We ig- 
nored data bins below 1.4 keV, where the MOS 
and PN spectra have significant discrepancies due 
to problems with the relative flux calibrations be- 
tween the two instruments. To better estimate the 
spectrum at low energies, we also fit the Chan- 
dra spectra of the central region from CXO 030720 
and CXO03082S (see the bottom panel of Figure 6). 


6 see the XMM- Newton calibration report CAL-TN- 001 8 - 2-3 


Because of the finer spatial resolution of Chandra , 
these spectra are mostly free from contamination 
by emission from the outer ejecta. Each XMM- 
Newton and Chandra spectrum was fit by a com- 
bination of a low and high temperature, optically 
thin thermal plasma model with independent ab- 
sorbing columns for each component, including a 
Gaussian line to account for the fluorescence iron 
line at 6.4 keV. The low temperature component 
accounts for the emission from the CCE. 

The parameters of the low temperature compo- 
nent were tied among all the spectra, and we as- 
sumed fixed solar abundances for all elements ex- 
cept for silicon and sulfur, whose abundances were 
allowed to vary simultaneously for all the spectra. 
We allowed the parameters of the hard component 
(which represents the variable emission of the stel- 
lar source) to vary among all the spectra, but we 
kept the abundances at the same (non-solar) val- 
ues in each spectrum. Table 3 shows the best-fit 
parameters for the constant component. The spec- 
trum can be fit with kT ~ 1.1 keV, Ah ~5.0 x 10 22 
cm“ 2 , with an absorption corrected luminosity of 
log Lx ~34 ergs s"" 1 , and low abundances for Si 
(~0.21 solar) and S (~0.47 solar). 

The model, however, was not accepted at above 
90% confidence for the following reasons: i) The 
Chandra spectra show significant excess below’ 
keV, perhaps originating in foreground emis- 
sion from the outer ejecta “bridge” (Weis et al. 
2004, see also Figure 2). The excess can be re- 
produced by including some flux from the outer 


11 




0.6 1 2 6 10 
channel energy (keV) 


Chandra Spectre during the X-ray Minimum 
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Fig. 6 . — Top: XMM-Newton EPIC 

PN spectra during the 2003 minimum 
(XMM 03 o722 ■ -XMMq 30818 ) after subtracting 

emission from the outer ejecta and Homunculus 
Nebula. Bottom : Chandra spectra of the central 
region in CXO030720 (black, HETG Oth-order) 
and CXO030828 (red, ACIS-S). The solid lines 
show the best-fit model of a simultaneous fitting 
of XMM-Newton (PN and MOS) and Chandra 
spectra by an individually absorbed 2T model. 
Barred, dot-bar, and dotted lines show the CCE 
component, variable component, and a Gaussian 
for the Fe K fluorescence line, respectively. 


ejecta using the model derived in §4.1.1. This 
slightly improves the reduced x 2 &om 1.51 for 
1205 to 1.50 for 1203 degrees of freedom. The 
included flux is ~0.4% of the whole outer ejecta 
emission, ii) The model overestimates the spec- 
tra near 2 and 2.7 keV, where lines of hydrogen- 
like Si and S ions are, and the model underes- 
timates the spectra near 1.8 and 2.5 keV, where 
lines of helium- like Si and S ions are. This suggests 
that the constant component includes plasma with 
temperature cooler than kT ~ 1 keV, though an 
inconsistency near 2 keV between XMM-Newton 
and Chandra spectra could also be caused by cali- 
bration uncertainties of the Chandra effective area 
around the iridium M-edge. Since cooler plasmas 
have lower emissivities in those lines, our simple 
IT model fitting for the CCE component can un- 
derestimate the elemental Si and S abundances, 
iii) The model underestimates the XMM-Newton 
spectra near 3—4 keV, possibly because the vari- 
able component is a multiple temperature plasma 
(see §4.2.2). iv) The Fe K line profile is complex 
and cannot be fit by a simple 1-temperature equi- 
librium model (see §4.2.3). Nevertheless, the pa- 
rameters given in Table 3 should approximate the 
physical properties of the CCE, except perhaps for 
the elemental abundances, which are more depen- 
dent on details of the model. 

As seen in Figure 6, the emission between 
1-3 keV in CXO 030720 comes predominantly from 
the CCE component. The 1-3 keV image in 
CXO 030720 was apparently point-like, which re- 
stricts the projected plasma size to within ~1 // 
of rj Car, the width of the Chandra point-spread 
function. This is equivalent to a projected distance 
<2300 AU at d ^2.3 kpc, suggesting that the hot 
gas which produces the CCE component is inside 
the Little Homunculus Nebula (which has a pro- 
jected extent of ~ ±2", Ishibashi et al. 2003). 

4.2* The Variable Component 

In order to determine the spectrum of the vari- 
able stellar component, we then subtracted from 
all the original spectra the recognized non- variable 
emission components - the outer ejecta, the X- 
ray Homunculus Nebula, and the CCE compo- 
nent. The remaining emission represents the vari- 
able component associated with the stellar source, 
and is shown in Figure 7. The top panel shows 
the variable component “corrected” for the con- 
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channel energy (keV) 

Fig. 7. — Top: XMM-Newton spectra of the variable emission. Bottom: XMM-Newton spectra renormalized 
to the line intensity of the helium-like Fe line. Best-fit models of the variable emission in XMMooor and 
CXO 030926 ? convolved with the EPIC PN effective area, are also displayed in grey and orange. Each color 
represents the same observation as in Figure 4. Boxes in the bottom panel show the high energy slopes of 
the bremsstrahlung emission at kT =4 keV with iVn = 0, 2, and 10 xl0 23 cm~" 2 (top), the bremsstrahlung 
emission at kT =1, 2, 4, and 8 keV without absorption (middle), and Fe line profile of the APEC model at 
kT =2, 4, and 8 keV with abundance of 0.5 solar (bottom). 
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stant components we identified. For comparison, 
we overlaid on the top panel of Figure 7 the best-fit 
models for the XMMooo? and CXO 030926 datasets. 

4.2.1 . Variation of the Overall Spectral Shape 

The bottom panel of Figure 7 shows the “cor- 
rected” variable components normalized at the 
Fe XXV line energy to highlight changes in the 
shape of the spectra. From XMM 0007 to XMM 0301 
(an interval of 2.5 years) the X-ray flux increased 
by a factor of 2-3 without remarkable changes 
in spectral shape. Four months later (observa- 
tion XMM 03 o 608 ) > the X-ray emission decreased 
only below 5 keV. Five days later (observation 
XMM 030613 )) when rj Car was close to the peak 
of the last flare before the 2003 minimum (Cor- 
coran 2005), the hard X-rays above ~3 keV in- 
creased by a factor of two. The X-ray emission 
dropped by almost two orders of magnitude by 
XMM030722, the first observation after the 2003 
minimum. In XMM 030802 ) the hard band flux 
recovered by a factor of 3 without any signifi- 
cant change in the soft band. The hard band 
flux stayed the same in XMM 030802 , XMMo30809 5 
XMMososis, and CXOo30828 (not shown), while 
the soft band emission slightly increased. The 
RXTE light curves (Corcoran 2005) also show a 
clear transition in the flux level and hardness ra- 
tio between XMM 030722 and XMMq 30802- This 
means that the X-ray minimum has two states. 
By CXO030926 5 only the hard band flux recovered 
to the pre- minimum level. 

The bottom panel in Figure 7 clearly shows 
that, throughout the observations, the hard band 
slope above 7 keV did not change significantly, nor 
did the ratio of hydrogen-like to helium-like Fe ion 
lines. This means that the electron temperature of 
the hottest plasma did not change during the 2003 
minimum. On the other hand, the soft band flux 
relative to the hard band decreased gradually from 
the X-ray maximum through the minimum to the 
recovery. This looks like an increase of the absorp- 
tion column to the X-ray plasma though the sit- 
uation is somewhat more complicated, as the fol- 
lowing sections describe. The 2003 minimum can 
be better described as an apparent decrease of the 
emission measure {EM.) as suggested by earlier 
ASCA observations (Corcoran et al. 2000), which 
means that either the amount of X-ray emitting 
material has declined, or that the amount that is 


visible to the observer is smaller. 

4.2.2. 1-Temperature Fit of the Entire Spectrum 

To quantitatively describe the spectral shapes 
and compare them with earlier results obtained by 
ASCA and BeppoSAX (see §5), we fit each spec- 
trum by an absorbed IT thermal APEC 7 model 
(which provided a better fit near the Fe K lines 
than the MEKAL model), with a Gaussian at 
6.4 keV to account for the Fe K fluorescence line. 
The results are listed in Table 4. 

Most of the fits were not acceptable at the 90% 
confidence level. One reason for the poor fit by 
the absorbed IT models is that the strong Fe Ka 
profile is quite complicated, as described in §4.2.3. 
Another reason, in particular for the spectra be- 
fore the 2003 minimum, is that the plasma is re- 
ally composed of multiple components with dif- 
ferent temperatures as shown by Corcoran et al. 
(2001b), as their best-fit plasma temperatures kT 
~4-5 keV do not account for the strong lines of 
helium-like S and Si ions seen in the XMM- Newton 
spectra (which should be emitted from plasma 
with kT ~ 1 keV). The plasma temperatures out- 
side the 2003 minimum [kT ~4-5 keV) are consis- 
tent with earlier results (Tsuboi et al. 1997; Cor- 
coran et al. 2000; Viotti et al. 2002; Leutenegger 
et al. 2003), but they are probably overestimated 
by our IT model fit (see the next section). The de- 
rived absorption column, which does not depend 
strongly on temperature for kT ~3—5 keV, grad- 
ually increased from V H ~ 5 x 10 22 cm“ 2 to 
^4 x 10 23 cm" 2 . 

4.2.3.. Fits to the E >5 keV Spectrum 

The EPIC spectra marginally resolved emission 
lines below ~5 keV and could not reliably mea- 
sure the neighboring continuum level. Spectral 
modelling of this energy band with unconstrained 
elemental abundances, therefore, has intrinsic un- 
certainty. On the other hand, the spectrum above 
5 keV is relatively simple and provides an unam- 
biguous measure of the hottest plasma from the 
continuum slope and Fe K line intensities. To 
characterize the hottest plasma, we thus fit the 
spectra above 5 keV by a simple absorbed IT 
model. 


7 http: / /cxc.harvard.edu/ atomdb/ 
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Table 4 

The best-fit models for the whole spectra 
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Fig, 8.— Spectra above 5 keV: EPIC PN (black), EPIC MOS1 (red) and MOS2 (blue), ACIS (green) 
and HETG (grey: -fist order, brown: -1st order). The solid lines show the best-fit absorbed optically 
thin-thermal plasma model (APEC code) with Gaussian components for Fe K a and p. 


















Fig. 8. — Continued. 


In each observation, we simultaneously fit all 
the available spectra, using an APEC thermal 
equilibrium model for the thermal emission. We 
simultaneously fit the spectra of two observations, 
CXO030720 and XMM 030722 J which were obtained 
during the 2003 minimum within 2 days of each 
other. A preliminary fit to the bright phase spec- 
tra showed that most spectra could be fit with 
kT =3.3 keV, a nickel abundance near 0.8 solar, 
with solar abundances for other elements. There- 
fore, we fixed the nickel abundance at 0.8 solar 
and the temperature at 3.3 keV, and we varied 
the Fe abundance, column density and emission 
measure. We included 2 narrow Gaussian lines to 
account for lines of Fe Ka and K/3 fluorescence, 
fixing their centroids at 6.4 keV and 7.1 keV, with 
their line widths (cr) fixed at 0.01 keV, and the K/3 
line flux fixed at 11.3% of the Ka line flux. We 
used photoelectric absorption cross-sections calcu- 
lated by Balucinska-Church & McCammon (1992). 
Table 5 gives the best-fit parameters, and the spec- 
tra, along with their best-fit models, are shown 
in Figure 8. The spectra except for XMM 0007 
and CXOXMM 0307 accept the assumed model at 


<90% confidence, mostly because these spectra 
show several spectral features which are not fit by 
the model. 

The most prominent of these features are 
excesses on the blue and red sides of the K- 
shell lines of helium-like iron. Those excesses 
are not caused by poor energy or gain cali- 
bration, since there is good consistency in pro- 
file between the XMM 030613 spectrum and the 
near- contemporaneous Chandra grating spectrum. 
The blue excess could be explained if the K 
lines of Fe XXV are Doppler broadened, with 
Av ~4000 km s’"" 1 . This is similar to the de- 
rived wind velocity of the companion star, v ~ 
3000 km s -1 (Pittard & Corcoran 2002). The red 
excess, especially during the 2003 minimum and 
CXOXMM0307 and XMM030802 when the X-ray 
flux was the weakest, would require a Fe XXV 
Doppler shift of Av ^7000 km s'" 1 , and such high 
velocities are not expected in rj Car. However, the 
red excess may be produced by unresolved emis- 
sion lines of Fe in ionization stages below Fe XXV 
in the 6. 5-6.6 keV band. Interestingly, the ASCA 
minimum spectrum in ASCA971224 showed a sim- 
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Table 5 

The best-fit models for the spectra above 5 keV 


Observation 

< fiorbit 

ZFe 

(solar) 

log 1 S.M. 
(cm" 3 ) 

EW(Fe) 

(eV) 

(10 22 cm- 2 ) 

reduced x 2 (d.o.f.) 

XMM 0007 

0.470-0.471 

0.48 (0.45-0.52) 

58.1 (58.1—58.2) 

87 (62-99) 

17 (14-20) 

1.43 (86) 

XMM0301 

0.922-0.924 

0.48 (0.47-0.50) 

58.5 (58.5-58.5) 

99 (85-104) 

17 (16-19) 

4.07 (82) 

XMMq 30608 

0.988 

0.55 (0.54-0.57) 

58.5 (58.5-58.5) 

164 (158-172) 

21 (20-22) 

4.21 (154) 

XMM03O613 

0.990 

0.44 (0.43-0.45) 

58.9 (58.9-58.9) 

119 (114-123) 

36 (35-37) 

3.66 (245) 

CXOXMMo 3 G 7 

1.008-1.009 

0.38 (0.30-0.47) 

57.0 (56.9-57.0) 

-233 (183-325) a 

40 (32-51) 

1.35 (90) 

XMM030802 

1.015 

0.50 (0.46-0.55) 

57.7 (57.6-57.7) 

<213* 

53 (48-58) 

1.91 (71) 

XMMo 30 S 09 

1.018 

0.70 (0.63-0.77) 

57.5 (57.4-57.5) 

<257 a 

32 (26-37) 

1.59 (54) 

XMM030818 

1.023 

0.62 (0.57-0.67) 

57.5 (57.5-57.5) 

<150 a 

27 (23-30) 

1.65 (84) 

CXOo 30828 b 

1.028 

0.33 (0.24-0.42) 

57.5 (57.4-57.6) 

-149 (96-214) a 

37 (28-47) 

1.81 (54) 

CXO 030926 

1.042 

0.33 (0.29-0.37) 

58.7 (58.6-58.8) 

110 (90-136) 

56 (50-62) 

1.54 (130) 


NOTE.— These fits refer to the spectrum after subtraction of the CCE and the Homunculus emission. The plasma temperature 
and the Ni abundance were fixed at 3.3 keV and 0.8 solar, respectively. 

a These values are regarded as upper-limits since the spectra did not clearly show the iron fluorescence line peak. 
b About 4% of the photon events pile-up, causing the spectrum slightly hard. 


ilar excess (see §5.2), while a Chandra grating 
spectrum near apastron (<f> ~0.53) showed a simi- 
lar broad- band excess (Pittard & Corcoran 2002). 
Presence of lower-ionization iron lines would indi- 
cate that the hottest plasma is no longer in colli- 
sional ionization equilibrium at 3.3 keV, and that 
the electron temperature is higher than the ion 
temperature. Non-equilibrium ionization (NEI) 
effects have been claimed in the X-ray spectrum 
of WR 140, another colliding wind system (Pollock 
et al. 2005). 

Similarly, a hump-like spectral feature is per- 
haps present between 5—6.4 keV and becomes 
more noticeable during the early phase of the min- 
ima. There are no specific emission lines in this 
energy band, except a weak Ca XIX line at 4.56 
keV. A similar feature is sometimes seen in the X- 
ray spectra of some AGNs (Tanaka et al. 1995), 
produced by gravitationally red-shifted material 
very near the central blackhole. However, there 
is little evidence that rj Car houses a black hole 
because of the lack of short term X-ray variability 
and relatively low X-ray luminosity. 

A second feature is a hard “tail” at >9 keV, 
best seen in the residuals in XMMooot> XMMo3oi, 
XMMo306Q8j and XMM 030613 (see Figure 8). This 
feature is probably not an instrumental or back- 
ground artifact: none of the spectra (except for 
CXOo 30828 ) suffer photon pile-up, and source 
count rates are much higher than the background. 


During the 2003 minimum the tail seems weaker 
or non-existent, though this may simply be an 
artifact since the source is weak at these times. 
The slope above 9 keV in XMM 030613 can be fit 
by a bremsstrahlung model with kT ~10 keV or 
power-law model with F ~ 2.3. Viotti et al. (2004) 
measured a similar photon index for a hard tail ex- 
tending up to 150 keV, seen in a BeppoSAX PDS 
spectrum in June 2000 (<t> ~0.46), but the flux 
between 13—20 keV, 1.4xl0” n ergs cm~ 2 s" 1 , is 
about three times larger than the extrapolation 
of the hard tail we measured in XMMo 306 i 3 - The 
hard excess we see in the XMM-Neivton data may 
be from an extremely hot plasma in the colliding 
wind region, or from non-thermal emission due 
to lst-order Fermi acceleration at the wind con- 
tact surface, which would produce a population of 
relativistic electrons that can upscatter UV pho- 
tons from the stellar photosphere and/or X-ray 
photons from the wind-wind shock. 

We also clearly see an excess at ~8.5 keV in 
the PN spectra in observations XMM 0301 and 
XMM 030613 * This feature might be due to emis- 
sion from K shell lines of heavy metals like Cu and 
Zn. Though the EPIC instrumental background 
shows emission from those elements (Striider et al. 
2001), inspection of the EPIC background data 
provided by the XMM project showed that the 
instrumental background at this energy was neg- 
ligible in spectra obtained outside the minimum. 
This feature could indicate an overabundance of 
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Cu or Zn, or perhaps be a part of the Cu edge at 
~9 keV (which would also require a large Cu over- 
abundance as well). On the other hand, the fea- 
ture might be produced by absorption edges from 
Fe XXV (~8.8 keV) and Fe XXVI at -9.3 keV 
(Lotz 1967). 

Finally, the fourth feature is a marginal dip 
around 7 keV, which is seen in all the PN spec- 
tra, but which is not so prominent in the MOS 
spectra. This dip could be caused by an underesti- 
mate of the line intensity of hydrogen-like Fe ions 
at 6.9 keV (perhaps due to an underestimate of 
the maximum plasma temperature) or an under- 
estimate of the strength of the Fe edge structure at 
7.1 keV, perhaps indicating larger Fe absorption. 
This dip is especially noticeable in the observa- 
tion CXOo30926? and will be discussed more fully 
in the analysis of our Chandra HETGS observa- 
tions (Corcoran et al. in preparation). 

4-2.4- Adding a Soft Component to Reproduce 
the Entire Spectrum 

Extrapolation of the model spectrum above 
5 keV from the XMMooo 7 j XMM 0301, XMM0306O8} 
XMM030613) CXO 030828 and CXO030926 data to 
lower energies underestimates the observed emis- 
sion at E <5 keV, indicating the presence of addi- 
tional cooler emitting material. We thus modeled 
the emission below 5 keV by fixing the best-fit 
model spectrum derived in §4.2.3, and adding an- 
other IT component with independent absorption. 
The results are given in Table 6 and Figure 9. 

These fits were not accepted at >90% confi- 
dence except the fit to the CXO 030926 spectra. 
One reason for the poor quality of the fits is that 
the fixed hard component, which reproduces emis- 
sion above 5 keV, somewhat overestimates the 
emission near 5 keV, which might indicate we 
slightly underestimated the absorption to the hard 
component. Another reason is that the best-fit 
models do not reproduce the strengths or loca- 
tions of many emission lines, in particular S XIV, 
which is generally emitted by plasma at kT > 1 
keV. This means that the spectra may require ad- 
ditional components with kT between 1—3 keV. 
The Ah values we derive from this 2 component 
modeling are very close to the ones from the IT 
fitting in §4.2.2. 


5. Comparison with the Previous Obser- 
vations Near X-ray Minima 

The weekly to daily monitoring of 77 Car with 
RXTE (Corcoran 2005) showed significant cycle- 
to-cycle variability in flux and X-ray . hardness. 
Since the RXTE data are contaminated by in- 
strumental and cosmic background, we compared 
our XMM-Newton and Chandra results with ear- 
lier spatially-resolved observations from ASCA 
(Tanaka et al. 1994), BeppoSAX (Scarsi 1993) 
and ROSAT (Triimper 1984). 

Between phase —0.10 < <f> 0 rbit < 0.10, there 
were 4 ASCA observations and 1 BeppoSAX ob- 
servation near the 1998 minimum, and 1 ROSAT 
observation during the 1992 minimum. These ob- 
servations are summarized in Table 7. Because 
of the extended point responses (FWHM ~a few 
arcminutes) of ASCA and BeppoSAX compared 
to XMM-Newton, source spectra extracted from 
these observations used relatively large (~ 3 ; ) ex- 
traction regions (Corcoran et al. 2000; Viotti et al. 
2002). This extraction area is about a factor of 
~20 larger than the area used in the XMM-Newton 
analysis (XMM SRC), and about a factor of ~5000 
larger than the area used in the Chandra analy- 
sis (CS SRC). Because the neighborhood of 77 Car 
is crowded with X-ray sources and diffuse emis- 
sion, the ASCA and BeppoSAX spectra will suffer 
greater contamination from neighboring sources. 

5.1. Comparison of the 1997/98 and 2003 
Spectra Outside the Minimum 

Outside the minimum, contamination of the 
stellar source spectra in the ASCA and BeppoSAX 
observations from the surrounding sources, the 
CCE component, and the reflection from the Ho- 
munculus Nebula is <10% above ~2.5 keV (see 
§5.2). 

Viotti et al. (2002) tried two models to fit their 
BeppoSAX spectra: an absorbed single temper- 
ature bremsstrahlung model -f Gaussian for the 
spectrum between 3—10 keV, and an absorbed 
two temperature model to fit the entire spectrum 
between M).l— 10 keV. Both gave consistent re- 
sults with each other for the central hard variable- 
component. On the other hand, in an analysis of 
ASCA spectra obtained just prior to the 1997/98 
minimum, Corcoran et al. (2000) fixed the emis- 
sion from the outer ejecta using the model derived 


19 




Fig. 9. — Fits of the spectra below 5 keV. 


Table 6 

Additional soft component for the spectrum below 5 keV 


Observation 

fiorbit 

kT 

(keV) 

Z 

(solar) 

log E.M. 
(cm~ 3 ) 

Nh 

(10 22 cm*” 2 ) 

Reduced x 2 (d.o.f.) 

XMM 0007 

0.470—0.471 

1.1 (1.1— 1.2) 

0.25 (0.20-0.31) 

58.3 (58.2-58.4) 

5.4 (5.2— 5.5) 

2.04 (225) 

XMM 0301 

0.922-0.924 

1.1 (l.l-l.l) 

0.18 (0.16-0.19) 

58.8 (58.8-58.8) 

5.4 (5.3— 5.4) 

8.85 (164) 

XMMo3060S 

0.988 

1.1 (l.l-l.l) 

0.24 (0.22-0.26) 

58.6 (58.6-58.7) 

6.9 (6.8— 7.0) 

5.41 (260) 

XMMo30613 

0.990 

1.6 (1.5— 1.6) 

0.25 (0.22-0.28) 

58.6 (58.5-58.6) 

7.7 (7.5— 7.8) 

12.02 (300) 

CXOo30S2S 

1.028 

0.5 (0.4— 0.9) 

<0.58 

59.9 (58.1-60.1) 

29 (17-39) 

1.32 (80) 

CXOq30926 

1.042 

0.9 (0.6— 1.8) 

0.28 (0.05-0.75) 

59.2 (58.0-60.5) 

23 (13-33) 

0.52 (109) 


Note. — These fits refer to the spectrum after subtraction of the CCE and the Homunculus emission. 
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Fig. 10. — Variation of the temperature and absorption column (Nr [whole]) vs. orbital phase. The filled 
marks are from the 2003 cycle (circle: XMM-Newton , triangle: Chandra , square: simultaneous fit between 
XMM-Newton and Chandra). The open marks are from the 1997/98 cycle (circle: ASCA , triangle: Bep- 
poSAX ). The ASCA spectrum during the 2003 minimum at phase 1.004 is consistent with the TVh and kT 
values from CXOXMMo 307 , as shown. 


Table 7 

Observations around the Minimum From Previous Cycles 


Abbreviation 

Satellite 

Date 

E orbit 

Exposure 

(ksec) 

N h 

(10 22 cm" 2 ) 

kT 

(keV) 

Referecnce 

ROSAT 920612 

ROSAT 

1992 Jun. 12, 22:33 

-0.998 

23.6 




A S C Ajj 70703 

ASCA 

1997 Jul. 3, 09:58 

-0.082 

12.6 

4.16±0.23 

3.79±0.31 

1 

ASCA 970719 

ASCA 

1997 Jul. 19, 15:11 

-0.074 

12.8 

4.09±G.30 

4.12±0.39 

1 

ASCA 971224 

ASCA 

1997 Dec. 24, 09:29 

0.004 

58.8 



1 

SAX980318 

BeppoSAX 

1998 Mar. 18, 05:01 

0.045 

39. 0 a 

15.4±0.4 

4.35±0.15 

2 

ASCA980716 

ASCA 

1998 Jul. 16, 04:40 

0.104 

13.4 

5.00±0.52 

4.89±0.S8 

1 


Note. — JVh and kT refer to the hottest components in 2- temperature fits of the entire spectra; ASCA results refer to 
analysis of the GIS2 spectra. Ref. l:Corcoran et al. (2000), 2:Viotti et al. (2002) 

a Exposure time for the MECS23 detector. 


in Corcoran et al. (1998). Though the coolest tem- 
perature they used to fit the ASCA spectra (kT 
~0.3 keV) is a bit lower than our adopted temper- 
ature (kT ~ 0.6 keV) and the resulting difference 
in the outer ejecta flux is about a factor of three 
at ~2 keV, this only produces a ~1% uncertainty 
in the flux of the variable emission outside the 
2003 minimum. The results of their BeppoSAX 
and ASCA analyses, therefore, are comparable to 
the result we derived from our IT fits to the entire 
XMM-Newton and Chandra spectra in §4.2.2. 

The plasma temperatures in the 1997/98 ob- 
servations were always around ~4— 5 keV, which 
is similar to our results. The derived values of Nr 


varied significantly with phase (Figure 10), reach- 
ing a maximum during the X-ray minimum, and 
declining thereafter. The variation in 2003 was 
similar to the Nr variation in 1997/98, except 
for the interval just at the end of the minimum. 
Though BeppoSAX and Chandra observed rj Car 
at a similar orbital phase (<£sax980318 ~ 0 , 045 , 
CXO030926 ^ 1.042), the Nr measured with Chan- 
dra in 2003 is a factor of two larger than the value 
derived by Viotti et al. (2002) from the BeppoSAX 
observations in 1998. This is consistent with the 
behavior of the hardness ratio after the recovery 
as measured by RXTE , which was higher in the 
2003 recovery than in the 1998 recovery (Corco- 
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Fig. 11. — Left : ASCA GIS2+3 spectrum on 1997 Dec. 24. Right : ROSAT PSPC spectrum on 1992 Jun 12. 
In each plot the spectrum of 77 Car during the X-ray minimum in 2003 is shown (red: the outer ejecta, blue: 
X-ray Homunculus Nebula, green: CCE component, purple: variable emission in CXOXMM0307; orange: 
surrounding sources within 3'. The solid black line shows the sum of these emission components. 


ran 2005). 

5.2. Comparison of the 1997/98 and 2003 
Spectra During the Minima 

In order to compare this dataset with the X- 
ray observations during the 2003 minimum, we re- 
analyzed the ASCA971224 GIS data. We screened 
the “revision2” data taken from the HEASARC 
archive 8 using gis clean and filtered it with the 
standard criteria, which excludes data affected by 
the South Atlantic Anomaly, Earth occultation 
and high background in regions with low geomag- 
netic rigidities. For both GIS2 and GIS3 detec- 
tors, we used a 3 ; radius circle centered on 77 Car 
as the source region SRC” , and used as back- 
ground emission from a 3' radius circle designated 
“3' BGD” centered to the north-east of 77 Car 
(Figure 2) where the soft diffuse X-ray emission 
from the Carina Nebula is apparently weak (Ham- 
aguchi et al. 2006). We used the standard GIS re- 
sponse files, version 4.0 (g[23]v4_0.rmf) and gener- 
ated ancillary response functions with ascaarfve r. 
3.10. We then merged both GIS2 and GIS 3 spec- 
tra together to improve signal to noise. The left 
panel of Figure 11 shows the GIS spectrum from 
ASCA 971224 * 

Because the stellar source was faint during the 
minimum and the ASCA point spread function 


8 http:/ /heasarc. gsfc.nasa.gov/W3Browse/ 


is large, contamination by nearby X-ray sources 
is significant for this observation. We therefore 
tried to estimate the contamination from emis- 
sion from sources within 3 7 of 77 Car “3 ; SRC” us- 
ing XMM030 72 2 (when the central point source of 
rj Car was the weakest) excluding the source re- 
gion “XMM SRC” used for the analysis of 77 Car 
from the XMM-Newton MOS2 data. We used the 
same background region as we used to estimate the 
ASCA background emission, “3' BGD”. About 
10% of the emission from 77 Car (which amounts 
to ~2Q% of the emission from the surrounding 
sources even in XMM030722) contaminates the “3' 
SRC” region. We simulated the contamination 
due to 77 Car and subtracted it from the ASCA 
source spectrum. The subtracted spectrum (Fig- 
ure 12) does not show any evidence of the ex- 
tremely strong nitrogen line at ~ 0.5 keV from the 
outer ejecta, suggesting that contamination is well 
removed. 

The resulting spectrum, after subtracting the 
simulated rj Car spectrum, shows weak emission 
lines at ~0.9 and 1 keV from Ne IX and Ne X, 
~0.8 keV possibly from Fe, and a marginal ex- 
cess around ~3~4 keV from Ca. This spectrum 
can be fit either by an optically thin, thermal 
(MEKAL) mpdel with 2-temperature components 
or by a cool component plus a power-law compo- 
nent for the hard emission at slightly below the 
90% confidence, assuming common absorption for 


22 





The spectral fit of the surrounding region 
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Fig. 12. — XMM-Newton MOS2 spectrum from 
XMM 030722 extracted from a 3' radius circle ex- 
cluding emission from r) Car. The solid line shows 
the best- fit 2T MEKAL model with common ab- 
sorption, while the broken lines show the individ- 
ual fit components. 

both components (Table 8). In either model, the 
average absorption to these sources is consistent 
w r ith the interstellar value 3 x 10 21 cm~ 2 , see 
§2.2 of Leutenegger et al. 2003). As the spectrum 
does not have strong emission lines, the elemental 
abundance of the 2T model is small (~0.04 so- 
lar) , while the abundances are not constrained by 
the IT plus power- law model. We thus fixed the 
abundances at 0.3 solar, a typical value for stel- 
lar X-ray emission obtained from low resolution 
CCD spectra (e.g. OB stars, Kitamoto & Mukai 
1996; Kitamoto et al. 2000; low-mass MSs, Tagli- 
aferri et al. 1997: low-mass PMSs, Yamauchi et al. 
1996; Kamata et al. 1997; Tsuboi et al. 1998.) 9 * * 
The observed flux is ^3.2 xlO" 12 ergs cm“ 2 s" 1 
(0.5-10 keV). 

The “3 / SRC 53 region includes at least 55 Chan- 
dra X-ray point sources (Evans et al. 2003): 
6 OB stars and 49 unidentified, weak sources 
that might be low mass pre-main-sequence stars. 
The combined absorption corrected flux from 
these sources calculated from Table 3 and 4 in 
Evans et al. (2003) is 1.3xl0“ 12 ergs cm" 2 s" 1 
(0.5—2.04 keV). There is also diffuse emission 


9 High resolution grating spectra with Chandra and XMM- 
Newton showed complex line profiles especially from OB 

stars, and therefore the abundance obtained from low res- 

olution spectra may not reflect real elemental abundance, 
(e.g., Cassinelii et al. 2001; Kahn et al. 2001; Miller et al. 

2002) 


within the source region which is also calcu- 
lated from Evans et al. (2003) to have a flux of 
^5. 7x10"“ 12 ergs cm" 2 s" 1 . The total absorp- 
tion corrected flux, ~7xl0'~ 12 ergs cm' 2 s" 1 , is 
roughly consistent with our result 8. 2 x 10" 12 ergs cm" 2 s' 
(0.5—2.04 keV, unabsorbed). Moreover, the com- 
bined absorption corrected fluxes of the point 
sources are as large as the flux from the hard 
component (1.7xl0" 12 ergs cm" 2 s” 1 ), and the 
flux from the diffuse emission is as large as the 
soft component, 6.5xl0" 12 ergs cm" 2 s" 1 , in the 
0.5 — 2.0 keV band. As seen in Figure 2, many 
point sources are detectable at intermediate en- 
ergies, where the hard component is dominant, 
but they are not clear in the soft band (except 
for HDE 303308). The hard component proba- 
bly represents emission from point sources with 
the soft component dominated by diffuse emission 
and HDE 303308. 

We have compared the ASCA spectrum to a 
model derived from our fits to the Chandra and 
XMM-Newton spectra including emission from 
the outer ejecta, from the Homunculus Neb- 
ula, from the CCE component and the variable 
model which fits the spectrum obtained from the 
CXOXMM 0307 observation, along with emission 
from the surrounding sources (diffuse emission and 
point sources) which fall within the 3 7 extraction 
circle. The ASCA spectrum agrees well with this 
combined model. Assuming that emission from 
the outer ejecta, X-ray Homunculus Nebula and 
the surrounding sources did not vary between 1997 
and 2003 and that spectral shapes of the CCE and 
variable components did not change dramatically, 
the flux of both CCE and variable components are 
the same within 50%. 

5.3. Comparison with the 1992 Minimum 

There is one ROSAT PSPC observation of 
77 Car during the minimum in mid-1992. Cor- 
coran et al. (1995) used this observation to show 
that the X-ray emission in the hard ROSAT band 
(E > 1.6 keV) decreased by a factor of 2 at that 
time. ROSAT however had very little effective 
area above 2 keV so that the bulk of the variable 
emission was not observable by ROSAT. 

Corcoran et al. (1995) extracted a spectrum 
from a 1/85 radius circle excluding a 30 // radius 
circle around the nearby bright star HDE 303308. 

To minimize contamination from surrounding 
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sources, we re-extracted the spectrum using an 
elliptical source region of 65" x 47." 5, which 
includes emission from the outer ejecta but 
which excludes HDE 303308 and point and 
diffuse sources around rj Car. We took the 
same background region as used in our anal- 
ysis of ASC A 971224 ♦ In analyzing the PSPC 

spectrum, we used the standard response file 
pspcb_gain2-256.rmf from the ROSAT calibra- 
tion databbase, but we generated the ancillary 
response file using the pcarf v. 2.1.3 as appropri- 
ate for this observation. 

The PSPC spectrum is shown in the right panel 
of Figure 11. We have overlaid all the spectral 
components of 17 Car on this spectrum as we did 
for the ASCA 971224 spectrum, except for the com- 
ponent due to the surrounding sources, whose con- 
tribution would be small in the ROSAT spectrum 
due to the smaller PSPC extraction region. Emis- 
sion between 0.5—1 keV, where the outer ejecta 
component is dominant, is about a factor of 30% 
larger than the model. This difference may be 
due to calibration uncertainties between ROSAT 
and XMM-Newton , or it might be produced by 
another variable source within the ROSAT extrac- 
tion area. Though the XMM-Newton model of the 
CCE component only contributes to the few high- 
est channels of the PSPC detector, this component 
is consistent with the ROSAT spectrum. Thus it 
appears that ROSAT did not detect the stellar 
component at all during the 1992 minimum. 

6. Discussion 

6.1. Emission from the Variable Compo- 
nent 

RXTE clearly showed that the variable emis- 
sion is periodic and leaves little doubt that the 
emission is produced by the WWC in a binary 
system (Corcoran 2005), and if so, that the X- 
ray minimum occurs near periastron passage. Our 
results further strengthen this conclusion. The 
temperature of the hottest plasma, ~3 keV, was 
mostly unchanged through the cycle, consistent 
with the collision of stellar winds at terminal ve- 
locities (e.g., see Ishibashi et al. 1999). Further- 
more, the spectral shape varied in a periodic way, 
as well. However, the 7% measured just after the 
end of the X-ray minimum was about a factor of 
two larger during the 2003 cycle than it was during 


the 1998 cycle. 

Most of the light curves derived from the Chan- 
dra and XMM-Newton data showed no short-term 
variations, suggesting that the X-ray plasma is 
produced by a steady source and is probably larger 
than Vshock x t exp , where Vs hock Is the pre-shock 
velocity and t exp is the exposure time of the obser- 
vation. For Vshock ~ 3000 km s -1 (corresponding 
to shock temperatures of ^ 3 — 4 keV) and expo- 
sure times of ^ 10 ksec (corresponding roughly to 
the XMM-Newton observing time) , the size of the 
emitting region is probably > 0.2 AU, and per- 
haps considerably larger. The RXTE lightcurve 
also shows clear variations in observations sepa- 
rated by as little as 1 day, suggesting that the size 
of the emitting region is 0.2 < r < 2 AU. Although 
the orbital elements are currently uncertain, the 
stellar separation at periastron may be as little as 
1.5 AU (Corcoran et al. 2001a), so that the size 
of the emitting region may be comparable to the 
stellar separation at periastron passage. 

Spectra outside the 2003 minimum require at 
least 2 temperatures, kT ^3 keV and 1-1.5 keV, 
consistent with the analysis of a Chandra grating 
spectrum by Corcoran et al. (2001b). The fail- 
ure to fit most of the spectra by 2 temperature 
models in §4.2.4 is perhaps an indication of the 
presence of more than 2 temperatures. This is 
consistent with a hydrodynamic simulation of the 
WWC plasma (Pittard & Stevens 1997), in which 
the plasma temperature is the highest near the 
stagnation point where the winds collide head-on, 
and cooler farther along the bow-shock. The ex- 
cess at the blue side of the Fe XXV line discussed 
in §4.2.3 may represent streams from the stagna- 
tion point flowing outward near the companion’s 
wind terminal velocity. 

The RXTE light curves also displayed quasi- 
periodic X-ray “flares” , which correspond to hard- 
ness ratio maxima (Corcoran et al. 1997; Corco- 
ran 2005). XMM 0306 O 8 occurred near the bottom 
of one of these flares, while XMM 030613 occurred 
near a flare peak. A comparison of these spectra 
in Figure 7 shows that the hard band flux went 
up by a factor of two, while the soft band flux be- 
low ^3 keV did not change except that emission 
lines are apparently weaker in XMMo 3 oei 3 - This 
is consistent with an increase of the hardness ra- 
tio measured with RXTE . This indicates either a 
slight increase in kT and/or E.M. for the hottest 
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plasma, possibly caused by increase of the density 
of the high-temperature gas near the stagnation 
point, as discussed for a similar event in the pre- 
vious cycle observed with A SC A (Corcoran et al. 
2000), The emission at E < 3 keV, on the other 
hand, did not vary strongly between those obser- 
vations. One possibility is that hard X-ray emis- 
sion from the highest-temperature plasma at the 
stagnation point of the shock cone varied strongly 
due to density fluctuations of the primary or sec- 
ondary wind, with the soft X-ray emission coming 
from a larger area which averages out such fluc- 
tuations. Another possibility is that the flares are 
related to instabilities in the WWC plasma. 

6.2. The Absorber around the WWC 

Plasma 

Before the 2003 minimum, Ah measured from 
our fits to the spectra above 5 keV (hereafter 
Ah[>5 keV]) was always about a factor of four 
larger than Ah measured from the spectra be- 
low 5 keV (Ah[<5 keV]; Figure 13). Ah[>5 keV] 
mainly represents Fe absorption, while Ah[<5 keV] 
represents absorption by lighter elements (Morri- 
son & McCammon 1983). Though this discrep- 
ancy can be produced by an Fe overabundance 
in the rj Car ejecta, the elemental abundances 
around the stellar source are near solar except for 
He and N, which are overabundant, and O and 
C, which are strongly depleted (Hillier et al. 2001; 
Verner et al. 2005; Davidson et al. 1984). Abun- 
dance anomalies of He, N, O and C would not 
change the amount of absorption above ^1 keV. 
This means that the soft emission does not suffer 
as strong absorption as the hard emission before 
the minimum, and hence this supports the idea 
that the variable emission originates in multiple- 
temperature plasmas with various absorptions. 
This is consistent with the idea that the hot emis- 
sion originates from a highly absorbed region near 
the apex of the shock cone, with the lower-energy 
emission arising from a region farther along the 
shock interface which suffers less absorption. Nei- 
ther Ah[>5 keV] nor Ah[<5 keV] varied signif- 
icantly between XMMooo? at 0 ~0.47 (close to 
apastron) and XMM0301 at 0 ~0.92, when the 
stellar separation changed by about a factor of 
3. The value of Ah[<5 keV] before the mini- 
mum (^5xl0 22 cm” 2 ) may be absorption beyond 
the wind of the primary. The higher value of 


Ah[> 5 keV] suggests that extra absorbing mate- 
rial must be located near the hottest plasma. 

After XMM0301 and through the onset of the 
2003 minimum, Ah[>5 keV] and the absorption 
column derived from fitting the spectra over the 
entire 1 — 10 keV band (Ah [whole], see Figure 10) 
increased to 3— 5xl0 23 cm“ 2 . As shown in Fig- 
ure 10 the observed maximum column density oc- 
curred near the mid-point of the X-ray minimum, 
and the column densities declined thereafter to- 
wards the end of the minimum. During the X-ray 
minimum, Ah[<5 keV] was unmeasurable because 
of the near total suppression of the low energy flux. 
This could mean that the E.M. of the cool compo- 
nent decreased during the minimum, or that the 
absorption to the cool plasma also increased on 
the line of sight. Nevertheless, the Ah[>5 keV] 
confirms that the X-ray minimum is closely asso- 
ciated with an increase in the amount of absorp- 
tion along the line of sight. This suggests that the 
WWC plasma entered into, or was hidden behind, 
the densest part of wind from the primary star. 
The behavior of the absorption towards the hard 
component, Ah[>5 keV], is in interesting contrast 
to the behavior of Ah [whole]. Ah[>5 keV] had al- 
ready increased to ~4xl0 23 cm” 2 by XMM 030613 
(i.e. prior to the start of the X-ray minimum) and 
it did not change strongly after that. This could 
be explained if the WWC region was bent, perhaps 
by the Coriolis force, causing the stagnation point 
to enter into the densest part of the primary’s wind 
earlier than the downstream material. 

The column density did not decline until some 
time after the 2003 minimum ended. This con- 
firms that the column density variation is not 
completely synchronous with the X-ray minimum. 
Conjunction probably occurred during the interval 
when Ah was near its maximum observed value, 
in the phase interval 0.99 < 0 < 1.05. We note 
that the Ah measured in 2003 (0=1.042) is a fac- 
tor of two larger than the Ah measured in 1998 
(0=0.045). Davidson et al. (2005, see their Fig- 
ure 2c) also noted that the absorption component 
in the P-?7(A9017) line in 2003.72 (0 « 1.0) was 
stronger than in 1998.21 (0 « 0.0), which might 
also indicate that the amount of absorbing mate- 
rial in the line of sight near the start of the X-ray 
minimum is increasing with time. However they 
also note that the equivalent widths of the Hr a 
emission lines were lower just after the end of the 
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Fig. 13— Variation of 2 Vh, EM. r intensity and EWof the Fe fluorescence line and RXTE fiux in - the 
2003 cycle. The XMM-Newton results are shown as circles. Chandra results as triangles r while squares show 
results derived from a simultaneous fit between XMM~Nevrton and Chandra. The dash line in the top panel 
represents our estimate of the absorption beyond the emit ral source. 




2003 minimum. This might suggest a distribution 
of wind material which has preferentially grown 
denser in the line of sight than perpendicular to \.rv. 
the line of sight in the 1998-^2003 interval. 

The equivalent width (EW) of the. Fe K fluo- . ' 
res.cence line varied between 100—200 eV T which 
is far smaller than the EW (~1.5 keV) produced 
by scattering X-rays from the Homunculus Neb- - 
ula (Corcoran et al. 2004). This suggests that the . 
X-ray emission , in all the observations comes di- 
rectly from the. WWC plasma, even during the 
X-ray minimum, with only a small scattered com- 
ponent (see Figure 14), From Inoue (1985), The ! 

Fe K fluorescence line EW produced by scattering' 
from a spherically-symmetric medium around an 
X-ray source is EW~75x A H [>5 keVJ/lO^eV,: ; ' 
where we have multiplied Inoue’s. coefficient by 
0.75, the flux ratio in the 7.1 — 9 keV band be- y 
tween a kT ^3.3 keV thermal emission component *• . 
appropriate to t he rj C ar spectra and the F = 1 . 1 . - . 

power-law emission component assumed by In— '- * 
oue (1985). In Figure 13, we plot the expected, y ^ •• 
Fe-K fluorescence line ■ ■■ EWs as open symbols-, for : 
a spiierically-symmetric -scattering medium . • The 7 
differences - between the . expected aiid observed *.v-- ' 
equivalent widths seem: to be larger when the . X- • 
ray flux goes up and smaller when it goes down 
(XM3 VIo30608-j XMM 030613 2 nd CXGo3G92s)> which ' 
may suggest -that the variation of the fluorescence % F 
line ’ is delayed by light Travel ■ time, b- . ' 

The- observed EWs are~25% smaller. than : : 
the expected- - EWs even: during The - • XMMooot •'-* . 
(<j*==QA70) and.'XIviMo 3 oi (<£=0.923) i observations;. ■ 
when- - the ■ direct X-ray ■ : emission ; did - r not .vary :• - ' ■- 
strongly.. This suggests That the -scattering, ma-. - 
terial is not spherically symmetric as assumed in’ 

Inoue (1985), and that there is a deficit of scat- 
tering material off the line of sight. This deficit 
is probably produced by the wind-wind collision,* 
since the region on the companion side of the 
shock has a much lower density than on the pri- 
mary star’s side. In order to reduce the fluorescent 
Fe emission by 25%, the full opening angle of the 
wind shock cone would need to be ^120°. This : 
is consistent with the measurement of the half- • 
opening angle of 50 — 120° by Pittard & Corcoran 
( 2002 ). 


6.3. The Nature of the X-ray Minimum : 

The 2003 X-ray minimum • can . be described b : 
as an apparent decrease in . E.M. r . along, with an > 
increase in N&. However, a column density of. 

Ah ~5x 10 24 cm' 2 is needed to reduce the emis- 
sion at 10 keV by two orders of magnitudes, 
and the observed maximum column density- (IV h 

5 x 10 23 cm"* 2 ) is much lower than this. The 
ASCA observations during the 1997/98 minimum. 
(Corcoran et al. 2000) also appeared consistent . 
with an apparent decline in E.M- 

Two mechanisms could produce the observed 
decrease in E,M. One possibility is an “eclipse 
. model 77 , namely that emitting region is almost to- 
t ally, but not completely, obscured by an • optically 
thick absorber with Ah> 10 25 cm" 2 , and we only' 
see the least absorbed residual emission during the 
minimum. Ishihashi et al. (1999) and Corcoran 
et aL (2001a) tested The eclipse model using the 
RXTE light curve of r) Car aiid showed that the 
duration of the X-ray minimum (~ 3 months) at 
high energies (JET > 5 keV) was difficult to explain 
for a spherically symmetric wind from : tj Car if the 
mass loss rate was constant at ~ 10~ 4 M© yr~~ l » 

The eclipse model seems to require that the wind is 
not spherically symmetric, or That the absorption 
to the hard X-ray emitting region: was enhanced 
for an extended interval near periastron passage. 

It also suggests that a portion of the hard emission b • 
is visible during the minimum so that the size, of T ?. 
the hottest region is greater than the size of the : 
occulting region. 

Another possibility is that the emissivity of the 
WWC itself fades during The minimum, i.e. that 
there is a decline in the amount of material hot 
enough to generate ~ 5 keV X-rays. Some WWC 
models predict a change in the X-ray emissivity : - 
of the shock around periastFon. For example, “ra- 
diative braking" (Gayley et al. 1997) in which UY 
photons from Car decelerate the companion’s 
wind at phases near periastron could reduce the 
intrinsic X-ray emission from the wind-wind colli- 
sion at E> 5 keV. Davidson (2002) suggested that 
strong instabilities near periastron might cause the 
shocked gas to radiatively cool, reducing the emis- . 
sion at high energies. Soker (2005) and Akashi : 
et al: (2006) suggested that Bondi-Hoyle accre- 
tion of the primary star’s wind by the companion 
star near periastron could suppress the compan- 
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Fig. . 14 — Cartoon of the 7? Car system; Winds from the primary star and the companion collide together; 
to form the wind-wind collision (WWC) region which has a conical shape. Behind the WWC boundary is a 
low density region filled with the wind of the companion star. The apex where the winds collide head on has 
hotter gas than the downstream material and also accumulates a buildup of cool primary wind material in 
front. We view the X-ray emission from the WWC plasma through the cool material at the apex and/or the 
undisturbed primary wind, and the absorber which lies far outside the binary system. Iron K fluorescence 
emission is produced by scattering iii the primary wind. One possibility for the CCE, where the CCE is 
produced in outer shock r is shown. Emission from the CCE only passes through the outer absorber beyond 
the binary systems - ‘ . / . - 


ion 7 s wind and thus reduce the emission from the 
colliding wind shock for a brief period: In opti- 
cal observationsy the He II A468F emission line, 
which is believed to arise near the WWC shocky 
decreased in concert with the X-ray drop (Steiner ' r 
& Damineli 2004)/ .. Martin et aL (2006) argued 04 .• 
that the disappearance of this line in a direct view 
of the star;, and in a reflected polar view (the so- 
called -F©S4vO position; Stahl et al. 2005) indi- 
cated a waning of the shocked gas during the min- 
imum. - 

Below, we examine some of the behaviors of the 
electron temperature; iron line profile, iron fluores- 
cence line variation and lack of variation at E > 

T keV during the minimum from the 2003 X-ray . 
spectra to try to decide between these two mech- 
anisms. - 

6.3.1. Electron Temperature 

The electron temperature of the hottest mate- • 
rial did not change significantly during the 2003 
minimum; as the normalized spectra at the bot- 
tom of Figure T~ clearly show. The constancy of 
the electron temperature of the hottest material 
argues against radiative braking, shock instabil- 


ity or accretion models which predict that the 
plasma temperature should decrease. Akashi et al. 
(2006) 7 s model in particular predicts complete 
shutdown of. the WWC activity during periastron 
passage, and therefore could only be consistent 
with a constant electron temperature if an alter-.: 
. native X-ray emission mechanism during the min- 
imum could produce plasma at nearly the same 
, temperature as the WWC emission. If the shock 
actually disappears during the minimum,, then the 
constancy of the high-energy X-ray continuum re- 
quires that a small part of the 3000 km s _1 sec- 
ondary wind continues to produce a shock by en~; 
countering material directly at normal incidence. 

6.3.2. The Iron Line Profile 

The variable excess emission on the red side 
of the Fe XXV" line can be explained by non- 
equilibrium ionization .(NEI) effects in . the high- 
temperature plasma, and variations in the strength 
of this feature means that the ion temperature 
of some of the high-temperature plasma in the 
colliding wind shock changed during the mini- 
mum. This might favor mechanisms in which the 
intrinsic emission from the shock declines with 
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time near periastron r though the NEI effect itself 
does not reduce emissivity of the continuum. The 
observed excess at ~6.55 keV is' consistent with 
emission lines of around Fe XXL present at an 
ionization time scale t‘~ nt ~ lO 10 ^ 11 s cm -3 • 
in the ~3.3 keV plasma,. where n is the gas den- 
sity and t is the time after the gas starts heating. 

Since the gas density n oc it should be larger. * 
by more than two orders of magnitude at perias- . 
tron than at apastron (using the orbital elements - ; 
of Corcoran et aL 2001a). Thus near periastron 
the plasma should reach ionization equilibrium 
much faster (~ 15 min, Soker 2005) than at apas- * • 
tron. Perhaps near periastron, both the cooling' 
time-scale and the equilibrium time-scale decrease 
relative to their values at apastron, but the cooling ■ 
time-scale decreases faster. This means that near 
periastron the fraction of non-equilibrium gas to 
gas in equilibrium is larger than it is near apastron. ' - 

However, since the radiative cooling timerscale has - ~ 
the same dependence on gas density as the ioniza- . * * 
tion equilibrium time-scale (oc n^ 1 ), an additional 
cooling mechanism (such as adiabatic expansion 
or conduction) needs to. be effective near perias- 
tron in order to decrease the cooling, time-scale 
faster than, the equilibrium time-scale. ; 

6.3.3. Iron Fluorescence Line Variations 

As noted in the previous section, the EWs ■ 

the Fe K fluorescence line stayed below ~2O0 eV % . 
during the . 2003' minimimy and therefore the ob- 
served X-ray emission is direct emissk>n r with only 
a small scattered component. The Fe fluorescence 
line flux decreased by a factor of ~50 during the 
2003 minimum (Figure 13). This could be pro- 
duced by reducing the number of E >T.IkeV X- • 
ray photons from the source seen by the fluorescing 
region, or by reducing the number of photons from 
the fluorescing region which reach the observer. If . 
we assume that the fluorescing region extends far 
from the stars then during the 2003 minimum ei- 
ther the X-ray source is nearly entirely covered by 
obscuring material along the line of sight to the 
fluorescing region, or the intrinsic X-ray emissivity 
has declined so that fewer ionizing photons reach 
the fluorescing medium. This picture is consistent 
with a similar variation of He II A4687 emission be- 
tween the direct emission from the central WWC 
region and polar emission scattered from the SE 
lobe (Stahl et al. 2005; Martin et al. 2006). 


6.3.4- The Lack of Variation at B > 7 he V Dur- 
ing the Minimum 

Tire hard band X-ray emission above keV * . 

reached its lowest value after the 2003 minimum 
started, recovered by a factor of 3 after XMM 030722 
and did not vary for at least ~3 weeks from 
XMM 030802 to CXOo 3082 S - Since RXTE light 
curves showed a similar trend m the previous cy- • ~ ' 
cle, the time variation seems to be phase repeat- 
able. The observed X-ray flux after XMM 030722 K 
was ~5% of the flux , at XMM 030613 * It is very 
difficult for an eclipse of the X-ray emitting re- 
gion by an opaque disk (similar to that proposed 
. by Ishibashi et al. 1999, for example) , to keep the 
observed flux constant for so long a time, since 
the stellar motion would be so fast near perias- 
tron that the emitting region moves beyond the oc- 
culter on ~day timescales, causing a short eclipse 
and a rapid rise in X-ray brightness. Perhaps a r 
“leaky absorber 77 eclipse, in which the absorbing 
medium is very extended and composed of opti-te - 
cally thick clumps, could better describe the be- . 
havior of the observed flux during the minimum. 
Models in which the shock cools may not be con- 
sistent with the constancy of the B : > TkeV emis— 
sion either, unless some residual emission at high 
energies could be generated after the shock cools . 

6.4* -The Origin o£ the GCJIL Component ^ 

The CCE component did not vary within t he ip - 
two months during the 2003 minimum when it 
could be observed. It apparently had the same 
spectrum during the 1998 minimum (as shown by 
our analysis of ASCA97-L 224)? an d did not show any ; 
enhanced emission during the minimum in 1992. 

This suggests that the X-ray source is stable on 
time-scales of 5— 10 years'. The Chandra image : 
during the 2003 minimum restricts the size of the 
emitting plasma to < 1", constraining the maxi- 
mum size of the. CCE at <2000 ALL The column 
density to this component, Ah ~5x 10 22 cm“ 2 , 
is smaller than the column density to the vari- 
able emission component around the 2003 mini- 
mum, but each component shows nearly the same 
column density in January 2003 (before the 2003 
minimum) and near apastron. After correcting for . 
absorption, the X-ray luminosity of the CCE is 
~10 34 ergs s-* 1 . .This suggests that, a stable lu- 
minous source is present outside of the absorbing 
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material which lies in front of the colliding-wind 
source at periastron. * 

One possibility is thafr the CCE component o 
is produced within the wind of the. companion 
star and/or t] Car due to the inherent insta- 
bilities of their radiatively-driven winds. The 
X-ray properties of the CCE component - kT 
— 1 keV with little apparent X-ray variation - re- “ 
semble the X-ray properties of normal OB stars i - , 
(e.g. Berghofer et ah 1997: Corcoran et al: 1993). 

For OB stars log £;*/£boi~ —6 to —8 (Berghofer S. 
et aL 1997), so that the X-ray luminosity of the 
CCE corresponds to a bolometric luminosity Lboi • 

— 10 4O ~ 42 ergs s -1 . This is somewhat larger 
than the luminosity of the companion according 
to Verner et aL (2005) (£*>©1 — 4 x 10 39 ergs s“ l ), .. 
but near the luminosity of p Car (Lboi — 2 x 
10 40 ergs s^ 1 , Hillier et aL 2001). 

Alternatively, the CCE emission could be pro- . v v 
duced in ambient gas shocked by a strong outflow 
from the stellar winds (see Figure 14). The maxi- * ' 
mum size of the CCE (r < 2000AU) means that it- - 
is inside the Little Homunculus nebula (Ishibashi • 
et al. 2003)^ so that the CCE component might be 
produced by the collision, of the polar wind from 
rj Car with the Little. Homunculus. The veloc- 
ity of the polar wind from r\ * Car is believed to 
be v — 1000 km s'” 1 (Smith et al. 2003a), which 
is sufficient to heat gas to kT ~ 1 keV. If all the 
kinetic energy of this outflow is converted to ther- h 
mal energy, the required mass loss rate would have 
to be only —3 xl0~ s Af e yr~* in order to match 
the observed emission measure (—710 X 10 56 cm" 3 );, 
which is only 3 x 10~~ 4 - 3 x 10~ 5 of the mass loss- 
rate of rj .Car (— 10^ 3 — 10~ 4 M 0 yrr\ Hillier et al. 

2001; Pittard & Corcoran 2002). The gas density 
inside the shell structure of the Little Homunculus 
should be less than the density of the shell, which 
is 10 6 < n h < 10 7 cm” 3 (Ishibashi etaL 2003). As- 
suming that radiative cooling dominates, the cool- : 
mg timescale of the plasma is r ~3nkT/n 2 A(T) r 
where A(T) is the emissivity. With a plasma den- 
sity of n <10 7 cm“ 3 . and plasma temperature kT 
— 1 keV, and A(lkeV) — 3.1xl0“ 23 ergs s'" 1 cm 3 ,: 
then r >0.5 years, consistent with our result. 

A third alternative is that the CCE might be 
associated with the outflowing material from the 
WWC that collides with cold circumstellar gas 
and is reheated. The required mass loss rate 
streaming out through the WWC is given by M == 


2Lx / eV 2 — 1 0 — / (c X (V/ 2000km s" 4 ) 2 ) M e yr^, 
where e is the efficiency to convert the kinetic en- 
ergy to X^ray enussion^and is the velocity of 
the outflow; which should be >1,000 km s” 1 ^ from 
the CCE temperatxire. Since the most of the com- 
panion’s wind goes into the shock cone, the M of 
the outflowing material should be — 10“ 5 M 0 yr" 1 
(Pittard .& Corcoran 2002), and e should only be 

T'J Conclusions &: Summary > 

We have observed 17 Car at key phases around 
the X-ray minimum in 2003 using the XMM- 
Newton and Chandra satellites and studied the 
variation of the X-ray spectral shape. . We also 
compared the behavior of the X-ray emission dur- 
ing the. 2003 minimum with observations during 
earlier minima made by. ASCA y BeppoSAX and 
ROSA T The results are summarized as follows. 

1. During the minimum, the hard and medium 
band emission from the central point source 
decreased to 1—5% of the flux outside the 
minimum, but the emission did not disap- 
pears The X-ray flux, varied in an apj>ar-: 
ently stepwise way near the middle of the 
minimum, suggesting that the X-ray mini- 
mum has two states. The light curve within 
any single observation did not show signifi- - 
cant variability on timescales of a few ksec, 
supporting that the hard X-ray emission is : 
produced in a WWC. 

2. The spectrum of the variable emission com- 
ponent changed markedly throughout the 
X-ray minimum; The medium band emis- 
sion decreased gradually relative to the 
hard emission. Assuming an absorbed IT 
model for the entire spectrum, the absorp- 
tion to the central source increased from 
5xl0 22 cm~ 2 to 4X10 23 cm" 2 ; The hard 
band flux went up prior to the X-ray min- 
imum, and decreased abruptly during the 
2003 minimum, without changing, the spec- 
tral slope significantly, suggesting that the 
plasma electron temperature is constant 
around 3 keV. 

3. The line flux ratios between the helium-like 
ions and hydrogen-like ions of Si, S, and Fe 
seen in the spectra before the 2003 minimum 


31 



indicate a similar type of multi-temperature 
plasma as seen near apastron by Corcoran 

et aL (2001b), ;-.v = 

4. The values of Afe derived from the spectra 
above 5 keV (JVh[> 5 keV]) are larger than 
those values of Ah derived from the entire 
spectra ( Afa [whole] ) or from the spectra, be— 
low 5 keY (iV H [<5 keV]) J but show less vari- 
ation than Nh [whole] or - A/fi [<5 . keV] . This v . 
suggests that the hotter plasma has stronger 
absorption, i.e. emission from tj Car is multi- 
temperature^ and each emission component 

is associated with an independent absorption; 
component. This means that the hotter and 
cooler emission originate in spatially distinct . 
regions. The hardest emission, which is be- ; • 

lieved to arise from material near the stag- . 
nation point of the WWC region^ perhaps k «... 
shows the effects of a build-up of material as 
t he colliding wind shock plows through the 
primary wind: / . * . • 

5. The Fe XXV line profile showed a small bliie 
excess, and a strong red excess before the -v 
2003 minimum^ The blue excess may be 
due to Doppler shifting of the line centroid 
or line broadening, and -if so, the observed 
velocity is ~4000 km s _1 . The red excess 
may indicate that some of the plasma is no ^:. 
longer in ionization equilibrium. The line 
showed a stronger red excess during the 2003V ' - 
minimum than outside the. minimum^ which . 
we interpret as decrease in the cooling time 
scale of the plasma near periastron. Inter- 
estingly, we note that an outburst spectrum . 
of the dwarf nova SS Cyg taken with the 
Suzaku satellite showed a similar red excess 

on the Fe K line profile accompanied with a 
decrease in the hard band X-ray ftux(Ishida. 
et.aJL 2007). ’ 

6. We find evidence for a hard tail at E > 9 
keY, at about an order of magnitude larger 
than the background levels in spectra out- 
side the X-ray minimum. This either in- 
dicates the presence of very hot plasma in 
the system or non-thermal emission perhaps 
produced by first-order Fermi acceleration of 
particles which then Compton-upscatter UV 
seed photons from the stellar photospheres 
and / or X-ray photons from the shock. 


7. The values of N& and kT of the X-ray 
spectrum at corresponding phases in pre- 
vious cycles are very similar to those during’: 
the 2003 cycle, except during the recovery ^ 
phase; In particular, the ASCA spectrum 
during the 1997/98 minimum (ASGA9n224) . 

is almost identical to the spectrum from 
CXOXMMosor* The difference in ATh during 
the recovery, as measured by Chandra and V 
BeppoSAX could be produced by an increase 
in the density of the wind from rj Car. 

8. The EWs of the Fe fluorescence line at 
6.4 keV are always below the values ex—. y- 
pected from Ah [>5 keV] assuming a spher- 
ically symmetric absorber. This is probably 
because the absorber is not spherically sym- 
metric and is consistent with the lack of 
scattering matter in the lower-density wind V: 
of the companion if the full opening angle of V 
the shock cone is ~ 120° T consistent with the * 
value derived by Pittard &r Corcoran (2002) 

9. The X^ray xniniminn is caused by the appar- * 

ent decrease of the plasma E~M~ An eclipse ; • 

model seems to explain the observed con-. : 
stancy of the plasma electron temperature; 
while the small EWs of the Fe fluorescence 
line and the deformation of the Fe XXV Be- 
like' triplet might indicate a real decline in - 
the amount of hot gas in the colliding wind 
shock. It might also indicate that nearly all; W 
of the shock is hidden by an optically thick 
but slightly porous stellar wind. 

10. We discovered a constant emission (GCE) 
component very near the central X-ray 
source. The GCE component can be seen 
between 1—3 keV during the X-ray mini- : 
mum; at other times the emission from the 
colliding wind source overwhelms it . This 
component showed no variation on either 
short timescales within any observation or 
between observations for about two months 
during the 2003 minimum, and was appar- 
ently the same during the 1997 / 98 and 1992 . 
minima. The GCE spectrum exhibited emis- 
sion lines of Si XIII and S XIV, and can 
be fit by an absorbed IT model with :2Vfr* 
~5xl0 22 cra~ z , kT ~ 1.1 keV and absorp- 
tion corrected Xx ~10 34 ergs •$"**, Chandra 
data suggest that the emission is point-like, 
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restricting the size of the source <2000 AU 
at 2.3 kpc. This emission might be produced 
by collisionaily-heated shocks from the fast 
winds from 7) Gar or perhaps the fast moving 
outflow from the wind-wind collision with 
the ambient gas, or might be associated with 
shocked gas which is intrinsic to the wind of 
r) Car. 

11. X-ray emission within 3' from r\ Car can be 
reproduced by 2 emission components with v 
interstellar absorption: cool plasma with kT 
~0.25 keV (perhaps coming from the diffuse . 
emission and HDE 303308), and hot plasma ~ 
with 2.6 keV or power-law F ~ 2.5, per- 
haps coming from a cluster of young stars.-. . 
Correcting for these emission components is 
important in understanding the X-ray spec- - 
trum of 7} Car derived from lower spatial res- . 
olution X-ray observatories like A SC A and 
Suzaku .. 


12. Emission from the outer ejecta did not vary* 
significantly on timescales of . several years. 

The slope of the highest energy emission / 
associated with the spectrum of the outer 
ejecta can be reproduced with ^Tr^O.G keVp 
which is consistent with the highest temper- 
ature derived from analysis of the emission 
line spectrum (Leutenegger et al. 2003). 

Currently, there are few reasons to suspect that rs; ; 
the variable X-ray emission is not driven by shocks 
produced by. the WWC from binary stars. For the 
first time, the spatially-resolved observations of 
7] Car obtained by XMM-N eurton and Chandra has - 
revealed the details of the plasma and absorption 
components which are directly associated with the 
WWC plasma. While the earlier spectral analysis 
of Pittard &, Corcoran (2002) did not account for 
the different emission components identified in this : * 
work, their fit was largely determined by emission 
above 2 keV which we have shown is dominated 
by emission from the WWC at the phase of their 
analysis. 

Still the biggest mystery is the cause of the X- /■ 
ray minimum r which is difficult to fully explain by 
occult at ion by a smooth wind from the primary 
star. If the X-ray activity really switches off dur- 
ing the minimum, a fading of the X-ray emission 
can be imaged in principle using the scattered X- 


ray emission from the Homunculus Nebula. The 
intensity of the Fe K-shell fluorescence line just af- 
ter theonset of the 2003 minimum should change 
differently between the eclipsing model and the 
fading model, and perhaps detailed monitoring of 
this line by spatially resolved X-ray spectral obser- 
vations could decide the issue. On t he other hand, 
deformation of the Fe XXV line profiles is unusual 
in the context of stellar systems, though similar 
deformation has been noted in iron lines arising 
from regions of strong gravity near black holes. - 
This feature requires further study with high res- 
olution spectrometers such as an X-ray micro- 
calorimeter to produce sufficient spectral resolu- 
' tion at.TkeV. 
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